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ASP:  A  Statistical  Package  in  APL 

[Draft  appendix  to  the  forthcoming  book, 

"Statistical  Computing  with  APL"] 

Designing  a  statistical  package  —  a  set  of  programs  for  statistical 
analysis  of  data  —  poses  problems,  seme  of  them  common  to  all  computing, 
some  perhaps  peculiar  to  statistics. 

Features  of  computing  for  statistical  analysis,  not  common  to  all  com¬ 
puting,  are  that  many  divers*  operations  are  carried  out  on  different  sets 
of  data,  that  several  basic  operations  may  be  applied  to  any  particular  set 
of  data,  and  that  only  after  a  baric  operation  has  been  applied  and  the  out¬ 
put  examined  can  the  user  see  (usually)  whether  the  operation  was  appropriate 
and  satisfactory.  Statistical  analysis  is  a  process  of  trial  and  error,  not 
to  be  specified  simply  and  completely  in  advance.  A  statistical  package  is 
therefore  n  vturally  thought  of  as  a  kit  of  tools  rather  than  as  a  single 
machine.  Ingenuity  may  be  devoted  to  organizing  the  tools  into  one  or  a  few 
machines,  but  if,  as  with  APL\36o,  the  computing  is  conversational,  flexibility 
is  to  be  preferred  to  grandeur,  unless  indeed  the  user  n^eds  to  be  shielded 
from  programing. 

A  pro-.r'itten  package  of  programs  is  attractive  to  a  user  who  does  not 
consider  hineelf  expert  in  the  language  in  which  they  arc  written.  Programs 
in  APL\36o,  such  as  those  presented  here,  demand  some  knowledge  of  the  language 
and  especially  some  lcnowledge  of  the  system  commands  a~d  of  how  to  manage 
workspaces.  Given  modest  knowledge  of  that  kind,  the  programs  can  be  success¬ 
fully  used  by  beginners  in  computing  who  would  not  readily  have  composed 
equivalent  progra  -.0  the  use  Ives.  Ultimately,  however,  tho  only  good  reason 
for  using  APL  for  statistical  work  is  to  have  complete  1 reedom  to  do  what 
one  wishes,  these  programs  should  be  thought  of  as  suggestions  for  how  to 
proceed  cad  as  a  source  of  examples  of  programming  technique,  rather  than  as 
prescriptions  for  good  statistical  practice.  The  programs  are  there  to  be 
altered  and  adapted.  They  are  net  locked. 

Programs  cm  have  at  least  four  technical  virtues,  apart  from  their  con¬ 
tent  ar.d  purpose,  namely  (i)  clarity  to  the  reader,  (ii)  speed  in  execution 
(CRJ  time  and  co-n.oct  time),  (iii)  economy  in  apace  required  for  execution, 

(iv)  eccnc-y  :n  space  required  for  storage.  All  four  are  important,  though 
their  relative  ivportance  is  arguable.  Any  actual  program  is  at  best  a  com- 
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promise  --  one  virtue  may  be  traded  for  another.  Only  if  a  program  is  to  be 
used  often  or  kept  a  long  time  are  the  economies  worth  much  thought,  once  a 
few  rudimentary  principles  have  been  learned:  loops  should  usually  be  avoided 
if  possible,  numerical  data  should  not  be  stored  permanently  in  floating¬ 
point  form,  and  things  like  that. 

Programs  (functions)  and  data  sets  to  be  kept  some  time  or  seen  by  others 
should  be  documented.  The  very  least  is  to  insert  a  line  or  two  of  comment 
in  the  definition  of  a  function  (as  for  example  in  RNORMAL  shown  below  in 
display-sheet  5  --  anyone  who  knows  what  is  meant  by  random  normal  deviates 
should  have  no  difficulty  in  using  this  function  after  seeing  a  print- out  ox' 
the  definition).  My  usual  practice  is  to  provide  minimal  information  about 
a  data  set  under  a  name  (of  a  variable  or  function)  stored  in  the  same  work¬ 
space,  consisting  of  the  letters  '■iKAT  joined  to  the  name  of  the  data  set, 
and  information  about  a  function  under  a  name  beginning  with  the  letters  HOT  . 
Such  documentation  occupies  least  storage  space,  usually,  if  it  is  a  variable 
consisting  of  a  character  vector.  Examples  are  displayed  below.  My  practice 
is  to  conclude  a  HCW-variable  with  a  date,  the  date  of  the  most  recent  per¬ 
ceptible  change  made  in  it.  The  user  may  quickly  discover  the  date  of  a  HCW- 
variable  by  calling  for  the  last  13  characters,  and  only  display  the  whole 
if  the  date  is  later  than  chat  of  his  last  print-out. 

Such  documentation  docs  not  attempt  to  explain  statistical  principles. 

More  extended  information,  explanations  of  purpose  and  illustrations  of  use, 
can  be  issued  as  duplicated  reading  matter,  not  kept  in  the  APL  system. 

Some  precautions 

As  far  as  conveniently  possible,  functions  should  be  written  so  that 
they  will  not  break  down  (be  suspended)  in  execution,  but  rather  will  ter¬ 
minate  with  a  v/aming  message  if  something  is  wi*ong.  At  least  to  begin  the 
definition  of  a  function  by  testing  the  arguments  (if  there  are  any)  for  com¬ 
patibility  seems  to  be  good  practice.  This  not  only  takes  care  of  a  common 
cause  of  failure  in  execution,  but  greatly  assists  a  reader  who  is  not  fully 
informed  about  vhat  the  function  does.  Tf  the  user  is  puzzled  to  get  the 
message  'NO  GO. '  when  he  tries  to  execute  the  function,  he  can  examine  the 
first  two  or  three  lines  of  the  definition  to  see  what  requirements  have 
been  placed  on  the  arguments. 
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Consider  for  example  the  xtinction  FIT  shown  below  (display-sheet  5). 

As  HCWFIT  explains,  this  function  is  intended  to  subtract  from  a  given 
vector  (Y,  say)  its  linear  regression  on  another  given  vector  (X,  say).  It 
could  be  used  to  detrend  a  time  series  by  removing  its  linear  regression  on 
the  time  variable.  The  means,  the  regression  coefficient  and  its  conven¬ 
tionally  estimated  standard  error,  are  printed  out  for  information.  The  number 
of  readings  N  is  defined  in  line  [l ],  and  thereafter  the  desired  calculations 
and  output  are  specified  in  lines  [3],  [ul,  [5],  [8l,  [9 1,  Do],  in  a  manner 
that  should  not  prove  difficult  to  read.  Line  [1 1  is  primarily  concerned  with 
testing  that  the  arguments  are  indeed  vectors  of  equal  length  not  less  than 
3,  and  if  they  are  not  line  [2]  causes  a  message  to  be  displayed  and  execu¬ 
tion  to  stop.  Line  [f]  tests  whether  all  members  of  the  first  argument  (the 
'independent"  variable)  are  equal.  If  they  are,  the  vector  X  is  this  line 
consists  of  equal  numbers  almost  or  exactly  zero,  and  the  regression  calcula¬ 
tion  in  line  [8]  should  not  (perhaps  cannot)  be  executed,  and  so  execution 
of  the  function  stops  at  line  [7].  If  execution  stops  at  line  [2],  the 
explicit  result  Z  is  empty.  If  execution  stops  at  line  [7],  the  explicit 
result  is  just  the  second  argument  minus  its  mean.  Only  if  both  these  traps 
are  passed  does  the  explicit  result  become  the  residuals  from  the  fitted 
regression,  as  intended. 

This  f’inction  FIT  has  not  been  protected  fPora  every  conceivable  cause 
of  breakdown.  If  the  arguments  were  vectors  of  equal  length  not  less  than  3, 
but  if  either  or  both  were  vectors  of  characters  rather  than  of  numbers, 
execution  would  be  suspended  at  line  [3]  or  [*']  where  the  sum  reduction  is 
called  for.  There  ought  to  be  a  primitive  monadic  function  in  API,  that  would 
yield  the  type  of  its  argument,  possibly  (for  example)  0  for  characters, 

1  for  single  bits,  2  for  4-byte  integers,  3  for  8-byte  floating-point  numbers. 
As  things  are  now,  one  may  test  whether  an  array  X  is  numerical  with  the 
expression 

r«iTcT,x  . 

This  works  because  APL\36o  distinguishes  empty  numerical  vectors  from  er.pty 
character  vectors,  but  that  has  the  appearance  of  being  an  accident  of  the 
implementation,  not  obviously  required  by  any  general  syntactic  principle; 
in  other  contexts  APL\36o  ignores  the  type  of  empty  arrays.  To  supply  a 
character  argument  to  a  function  where  a  numerical  argument  is  intended,  or 


-4  - 

vice  versa,  is  no  doubt  an  unusual  error  in  execution,  scarcely  worth  guard¬ 
ing  against.  At  any  rate,  it  has  not  been  guarded  against  here. 

Some  of  the  functions  in  this  collection,  such  as  DCWNPLOT  ,  require 
that  one  or  more  global  variables  should  have  been  defined  before  the  function 
is  executed.  There  seems  to  be  no  way  at  present  to  test  in  the  function 
definition  whether  that  is  so,  without  a  suspension  if  in  fact  the  required 
global  variables  have  not  been  defined. 

Another  possible  reason  for  breakdown  in  execution  of  a  correctly  written 
function  is  that  there  is  insufficient  rocm  in  the  active  workspace.  If  the 
space  requirement  vrere  determined,  as  a  function  of  size  of  arguments  and 
any  other  relevant  variables,  this  also  could  be  tested  for.  I  have  generally 
not  done  so,  but  an  example  of  a  possibility  of  this  sort  may  be  seen  in  the 
function  FILTER  (display-sheet  9).  The  meat  of  the  function  is  in  the 
single  line  [9],  which  executes  fast  but  needs  room,  sometimes  a  lot  of  room. 
The  last  three  lines  in  the  definition  constitute  an  alternative  program  to 
lir>*»  01  slower  to  execute  because  of  the  loop  but  needing  less  space.  In 
line  [8]  a  test  is  made  of  available  space,  with  a  branch  to  [10]  if  [9] 
cannot  be  executed.  (For  this  calculation  it  is  supposed  that  the  argument 
X  consists  of  floating-point  numbers,  which  it  generally  does  in  practice.) 
Breakdown  in  execution  because  of  lack  of  space  has  not  been  guarded  against 
completely,  there  may  be  too  little  rocm  in  the  active  workspace  even  for 
the  mere  economical  program,  but  the  worst  space  problem  in  FILTER  has  been 
taken  care  of. 

Suspensions  of  function  execution  will  occur.  Sometimes  the  user's 
attention  is  so  fixed  on  remedying  the  causa  or  taking  other  action  prompted 
by  the  suspension  that  he  omits  to  terminate  execution  of  the  suspended  func¬ 
tion.  Thenthe  function's  local  variables  supersede  any  global  objects  with 
the  same  names.  The  most  annoying  effect  of  this  kind  comes  from  statement 
labels,  which  cannot  (during  suspended  execution  of  a  function)  be  treated 
like  ordinary  names  of  variables  and  assigned  new  values.  Preferably  state¬ 
ment  labels  should  be  unlike  any  other  names  used.  Accordingly  all  statement 
labels  here  begin  with  an  underscored  L  .  If  the  user  refrains  from  ever 
assigning  to  a  variable  or  function  a  name  beginning  with  underscored  L 
(no  very  irksone  constraint,  since  underscoring  is  troublesome  anyway), 
there  can  be  no  confusion. 
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Possibly  there  should  be  some  similar  distinguishing  of  all  local  vari¬ 
ables  in  defined  functions.  That,  however,  is  not  so  easy  to  achieve,  partly 
because  in  developing  a  function  the  author  should  feel  free  to  change  the 
status  of  variables  between  local  and  global.  For  economy's  sake,  no  local 
variable  should  have  a  name  longer  than  three  characters*  .  If  global  vari¬ 
ables  and  functions  intended  for  permanent  storage  are  given  descriptive 
names  more  than  three  characters  long  (as  is  ray  usual  but  not  invariable 
practice ),  at  least  there  will  be  no  confusion  between  them  and  local  vari¬ 
ables. 

A  particular  kind  of  suspension  occurs  with  quad  input,  or  as  the  User's 
Manual  calls  it,  evaluated  input.  Only  three  of  fhe  functions  shown  here  have 
this  feature:  INTEGRATE,  SCATTERPLOT,  TSCP.  In  answering  the  first  question 
in  SCATTERPLOT  or  TSCP  ,  the  user  will  normally  type  two  numbers,  as  sug¬ 
gested;  but  if  instead  he  gives  the  name  of  a  previously  defined  2 -element 
vector,  this  must  be  different  from  the  names  of  all  eleven  local  variables 
(not  counting  statement  labels)  in  the  function.  The  same  applies  to  the 
third  question  in  TSCP  ,  where  now  there  may  be  greater  temptation  to  reply 
with  the  name  of  a  previously  defined  character  vector  (Z  would  be  permis¬ 
sible,  or  any  name  with  more  than  one  character).  Similarly,  if  the  function 
asked  for  in  INTEGRATE  is  previously  defined,  its  name  must  not  be  any  one 
of  A,  H,  IT,  7  (it  could  happily  be  F,  or  any  name  with  more  than  one  character). 


'  There  is  a  small  unobvious  economy  in  storage  of  functions  to  be  derived 
from  keeping  the  names  of  local  variables  short,  mostly  single  characters.  That 
is  that  in  a  collection  of  functions  the  same  names  of  variables  will  occur 
again  and  again,  for  there  are  only  26  letters  in  the  alphabet  to  choose  from. 

In  storage  each  different  name  appears  just  once  in  the  symbol  table,  which 
therefore  needs  to  be  large  enough  to  accommodate  all  the  different  names  in 
the  workspace  —  the  names  of  groups,  functions  and  variables,  including  the 
variables  and  statement  labels  in  the  function  definitions.  In  the  functions 
shown  here,  a  handful  of  cr.e-letter  names  account  for  many  of  the  local  vari¬ 
ables;  each  appears  several  times  in  the  definitions,  but  onlj  once  in  the 
symbol  table.  Therefore  these  functions  can  be  stored  in  workspaces  with 
symbol  tables  substantially  reduced  below  the  standard  setting. 
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Miscellaneous  remarks 

Perhaps  a  word  should  be  said  about  some  of  the  algorithms,  with  reference 
to  their  speed  in  execution.  One  should  bevrare  of  supposing  that  comparisons 
of  alternative  algorithms,  executed  in  compiled  Fortran,  are  relevant  to 
apparently  similar  algorithms  executed  in  APL\36o.  The  richer  vocabulary  of 
primitive  functions  in  APL,  and  the  different  handling  of  arrays,  can  lead 
to  different  conclusions.  Because  of  its  brevity,  the  function  RNORMAL 
seems  to  be  about  the  fastest  possible  generator  of  random  normal  deviates 
in  APL.  The  function  HA  is  a  slow  Fourier  transform,  but  performs  well 
with  series  of  modest  length,  such  as  can  be  accommodated  in  the  customary 
32K  byte  workspace ;  and  it  leaves  the  user  completely  indifferent  to  the 
prime-power  factorization  of  the  iength  of  the  series.  (A  fast  Fourier  trans¬ 
form  will  no  doubt  eventually  be  added  to  the  collection.  )  The  loop  at  the 
end  of  SCATTERPLOT  could  be  avoided  by  an  outer  product,  but  that  teems 
slow  in  execution  and  certainly  requires  considerable  space  --  a  counter¬ 
example  to  the  general  principle  of  avoiding  loops  if  possible. 

Most  of  these  functions  have  been  developed  independently  of  anyone 
else's  programming.  An  exception  is  INTEGRATE  ,  based  on  a  program  for 
integration  by  Simpson's  Rule  due  to  John  Richardson  (APL  Quote  Quad,  2/3 
(1  970),  26-27  ).  The  function  JACOBI  was  developed  from  two  APL  programs 
written  by  others,  but  is  not  now  very  close  to  either  source.  K.  W.  Smillie's 
collection  of  statistical  programs,  STATFACK2  ,  provided  many  interesting 
examples  of  APL  technique,  but  I  have  deliberately  refrained  from  actually 
using  his  programs  in  order  not  to  be  influenced  in  statistical  thinking. 

The  display-sheets  following  show  function  definitions  and  documentary 
variables  in  the  workspaces  ASP  and  ASP2  .  In  alphabetical  order,  with 
the  sheet  number  on  which  they  appear,  they  are: 

Jo  'umentary  variables 


DESCRIBE 

HOU JACOBI 

'7 

JCWC  OUT  IITGENC  Y 

4 

HC  .’MAX 

?4 

HOWFILTER 

0 

HOI /PLOT 

13 

HOWFIT 

3 

HOi  (REGRESSION 

4 

HGWHARMONIC 

1  0 

110! /SUMMARIZE 

2 

HOUHUBER 

1  2 

HO' 'TRIPLE  PLOT 

22 

HCN'INTEGRATE 

2! 

tai’nifiittiytf  Mrmrfi^yiiirtiii'fniii^frr  iiiniMiiiiiirrf  n 
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Function  definitions 


ANALYZE 

8 

AUTOCOV 

1 1 

CONTINGENCY 

15 

DOWNPLOT 

20 

EFFECT 

8 

FILTER 

9 

FIT 

3 

FOURFOLD 

l  6 

HA 

1  0 

HUHER 

l  3 

HUBER1 

1  3 

HUBER 2 

1  3 

IN  IF 

1 1 

INTEGRATE 

21 

JACOBI 

17 

MAV 

1 1 

MAX 

2  4 

multipoly 

16 

NIF 

1 1 

^VWVAP*.*1H  ■  ■  -W *•*■*■»&*>. 


POOL 

REGR 

REGRINIT 

RNORMAL 

ROWCOL 

ROUCOLDISPLAY 

ROWCOLFERMUTE 

ROVJPLOT 

SAMPLE 

SCATTERPLOT 

SHOP' 

STDIZE 

STRSS 

SUMMARIZE 

TAPER 

TDP 

TSCP 

VARIANCE 
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w 

WfI 


)LOAD  ASP 

SAVED  21.48.43  07/17/73 
)FiJS 

ANALYZE  DOUSE  LOT  EFFECT  FIT  IMF  INTEGRATE  JACOBI  UIF 

REGR  REGRINIT  RNORMAL  ROUCOL  ROUCOLDICPLAY  ROUCOLPER14UTE  ROUPLOT 

SAMPLE  SCATTERPLOT  SHOW  STDIZE  STRES  SUMMARIZE  VARIANCE 

)VARS 

DESCRIBE  tlOWFIT  SOW  INTEGRATE  iiOUJACQDI  HOUPLOT  HOUREGRESSIOU 

HOWSUMMARI  ZE  A  IL  Q. 

)GRPS 

imYGROUP  REGRESS IOUGROUP  ROUCOLGROUP 

)GRP  NWAYGROUP 
AiJALYZE  EFFECT  £ 

)GRP  REGRESSIONGROUP 

REGR  REGR  IS  IT  SHOU  STRES  VARIANCE  Q 

)GRP  ROUCOLGROUP 

ROUCOL  ROUCOLDISPLAY  ROWCOLPERMUTE  £ 


DESCRIBE 

ASF  -  A  STATISTICAL  PACKAGE  IE  A  PROGRAMMES  LANGUAGE 

t++  .  t  t :  + . : : . . :  . .  * 

SONIC  PROGRAMS  FOR  STATISTICAL  ANALYSIS  BY  F.  J.  AUS COMBE,  WHO  WILL  APPRECIATE 
COMMENTS  AND  ERROR  REPORTS.  {PHONE:  203-432-4752.) 

' HOUPLOT '  AND  ' HOWTRIPLEPLOT '  DESCRIBE  PLOTTING  FUNCTIONS . 

SEE  ' HOWSUMMARI ZE'  FOR  SIMPLE  SUMMARY  STATISTICS,  ' HOWFIT '  AND  ' HOW REGRES¬ 
SION  '  FOR  LEAST-SQUARES  REGRESSION,  AND  ' HOUHUDER '  FOR  ROBUST  REGRESSION.  SEE 
ALSO  THE  FUNCTIONS  ' ROUCOL '  AND  'ANALYZE'.  ' HOWFILTER'  AND  ' HOUHARMONIC'  DES¬ 
CRIBE  FUNCTIONS  FOR  ANALYSIS  OF  TIME  SERIES ;  SEE  ALSO  'AUTOCOV' .  ' HOUCONTINGEN- 

67'  DESCRIBES  FUNCTIONS  FOR  ANALYSIS  OF  CONTINGENCY  TABLES. 

' NIF '  AND  ' IN  IF '  ARE  THE  NORMAL  INTEGRAL  FUNCTION  AND  ITS  INVERSE.  VARIOUS 
MATHEMATICAL  OPERATIONS  ARE  DESCRIBED  IN  ' UOWJACODI' ,  ' HOW INTEGRATE ' ,  'HOUMAX' . 

THE  FUNCTIONS  'ROUCOL',  'ANALYZE',  'EFFECT',  'AUTOCOV',  'SAlfPLL" ,  'RNORMAL', 
' NORMALIEST'  ARE  DESCRIBED  IN  CHAPTlR  7  OF  STATISTICAL  COMPUTING  WITH  APL. 

THESE  FUNCTIONS  RESIDE  IN  WORKSPACES  NAMED  'ASP'  AND  'ASP2'  .  SOME  OF  THE 
FUNCTIONS  ARE  ORGANIZED  IN  GROUPS.  BEGIN  NY  LISTING  ALL  FUNCTIONS,  VARIABLES 
AND  GROUPS  IN  EACH  WORKSPACE. 

1 -ORIGIN  INDEXING  IS  ALWAYS  SUPPOSED.  SOME  FUNCTIONS  NEED  THE  GLOBAL  SCALAR 
CHARACTER  VARIABLE  C  CAUSING  A  CARRIER  RETURN.  C  IS  INCLUDED  IN  EVERY  GROUP. 
C  SHOULD  BE  SPECIFIED  OR  COPIED  WHEN  ANY  OF  THESE  FUNCTIONS  {BUT  NOT  A  GROUP)  IS 
COPIED.  THE  SCALAR  VARIABLE  B  CAUSES  A  BACKSPACE.  THE  VECTOR  4  IS  THE  AL¬ 
PHABET,  FIRST  PLAIN  THEN  UNDERSCORED. 

DEPT.  OF  STATISTICS,  YALE  UNIV . ,  HE W  HAVEN,  CT.  06520  -  11  JULY  1973 
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ten***  *m+r..+.  .  . 


UOUSUMARIZE 

SUM-ARY  STATISTICS  OF  A  DATA  SET 
SUMARIZL  I 


THE  ARGUMENT  (Z)  MAY  HE  ANY  NUMERICAL  ARRAY  HAVING  AT  LEAST  4  MEMBERS.  MEAS¬ 
URES  OF  LOCATION,  SCALE  AND  SHAPE  OF  DISTRIBUTION  ARE  PRINTED  OUT.  SHAPE  IS  IN¬ 
DICATED  BY  FREQUENCIES  OF  THE  DATA  IN  THE  SIX  INTERVALS  INTO  WHICH  THE  REAL  LINE 
IS  DIVIDED  BY  THE  MEAL'  +  "21012  *  TilE  STANDARD  DEVIATION  OF  THE  DATA.  EX¬ 
PECTED  COUNTS  IN  THESE  INTERVALS  FOR  A  RANDOM  SAI-IPLE  OF  THE  SAME  SIZE  TAKEN  FROM 
A  NORMAL  DISTRIBUTION  WITH  THE  SAt-IE  MEAN  AND  VARIANCE  ARE  GIVEN  FOR  COMPARISON . 
SHAPE  IS  ALSO  INDICATED  BY  FISHER'S  SKEWNESS  AND  KURTOSIS  MEASURES,  WITH  THEIR 
STANDARD  ERRORS  FOR  A  RANDOM  SAMPLE  OF  THE  SAME  SIZE  FROM  A  NORIAL  DISTRIBUTION. 

28  JUNE  1973 


V  SUMMARIZE  Z-,A;D-,K;M;N;S-,N 

Cl  3  -*2+4</l/'«-pZ«-,Z 

[2]  -t-pO 'THE  ARGUMENT  MUST  BE  AN  ARRAY  WITH  AT  LEAST  4  MEMBERS.  ' 

[3]  C,' NUMBER  OF  READINGS:  '  \N 

[4]  'EXTREMES:  '  ;(Z«-ZUZ])[1 ,4] 

[5]  '  QUARTILES :  '  .  5*  ( +/Z[ft'+1-/1  ] ) ,+  /ZU<-(  L 1+A  *4),r</*4] 

[6]  'MEDIAN:  '  ;0.5x+/Z[  ( Ll+W+2)  J//+2] 

[7]  'MEAN:  '  NM+/Z) HI 

C  8  ]  Q , '  INTERQUARTILE  RANGE :  '\-/A 

[9]  'VARIANCE:  '  ;i(<-(S++/W*-Z*Z<-Z-l!)*N-l 

[10]  'STANDARD  DEVIATION:  ' 2 

Cll]  /K0.1*C0.5+0.0001x/1/x  2275  13591  34134 

Cl  2  D  FREQUENCIES:  '  ;+/(*4+ 16  )<> .  =2C~3r  LZ+C; ' ,  WITH  FITTED  NORMAL  EXPECT  AT  10 

NS  '  ;4,<M 

C13]  'SKEWNESS  MEASURE  Gl:  ';( iK*D)*  */{+/!  ’*Z )  ,/'+  "l  0  ~2  WITH  NORMAL  S.E. 

'  ;(*/6,//+  "2  0  1  "1  3 )*+  2 

C14]  iK*K ) x ( +/ ( +  /r/xf / ) , Nr  "1  0  "2  1  "3)-3xJx5t(/;-2)x//-3 

Cl 5]  'KURTOSIS  MEASURE  G2:  '  ,A  ; '  ,  WITH  NORMAL  S .E .  ' ;(  */24,A'+  _3  C 

"2  “1  3  "1  5)**2;£ 
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ilOWFIT 

SIMPLL'  REGRESSION 
Z*-X  FIT  y 

THE  TWO  ARGUMENTS  MUST  HE  VECTORS  OF  EQUAL  LENGTH  NOT  LESS  THAN  3. 

THE  MEANS  OF  THE  TWO  ARGUMENTS ,  AND  THE  REGRESSION  COEFFICIENT  OF  THE  SECOND 
ARGUMENT  ON  THE  FIRST,  WITH  ITS  CONVENTIONAL  ESTIMATED  STANDARD  ERROR,  ARE  DIS¬ 
PLAYED.  THE  EXPLICIT  RESULT  IS  THE  VECTOR  OF  RESIDUALS. 

30  MAY  1973 


V  l+X  FIT  Y\U\U\S,U\V 

tl]  -*-2+(tf=+/pY)A(3^^/py)A(l=ppy)Al=ppy 

C2]  -*pO«- 1  ARGUMENTS  MUST  BE  VECTORS  OF  EQUAL  LENGTH  NOT  LESS  THAU  3. '  ,Z«-' ' 

[3]  x+x-u+(+/x)*n 

[4]  Z<-Y-VH+/Y)iN 

[5]  'MEANS  ARE  *  ;£/.  V 

[6]  -*-8-A/X=lt;r 

[7]  -pQ-'/VO  REGRESSION  CALCULATED.' 

[8]  'REGRESSION  COEFFICIENT  IS  '  iB+-(X+.*Z)iS<-X+.*X 

[9]  5^(Z+.xZ^Z-flK^)fSx/V-2 

[10]  '  WITH  ESTIMATED  STANDARD  ERROR  '  ;S* 0.5 ;  •  (  '  ill- 2 ; '  D.F.)' 

V 


V  Z+RNORMAL  N‘,V 

[  1  ]  -*-2+ v  /  (/ite  1 )  a(/V= IN )  a  1  =p  N+ ,  11 

[2]  -*-0.pL>,/70  GO.  '  ,Z+-'  ' 

[3]  Z«-(*2147483647)*?(2,nV*2)p2147483647 

[4]  Z«tft(,  1  2  0  .oZ[  2 ;  ]xo2  )xK,  !/•♦•(  _2X#Z[1 ; ] )*0. 5 

[5]  ft  fl/WDOM  NORMAL  DEVIATES  NY  DOX -MULLER  METHOD. 

V 


V  Z+i/  SAI'JPLE  CP 

[1]  -^2+(  A/,  (/•/  =  [// )Al</;)A(l2pp/^)AA/(l>C,P)A(.)<CP^-,CP 

[2]  ->p[>V0  (7(7. ' 

[3]  Z*"+  ^CP°  .  <(  *2147483647  )  * ?Pp 2 14740  3647 

[4]  pi  SAMPLE  OF  SIZE  11  FROM  A  DISTRIBUTION  OVER  NOFNEGATIVE  INTEGERS  HAVING  C 
UMULATIVE  PROBABILITIES  CP. 


V 


NGWRCGRLSSLOt! 


MULTIPLE  REGRESSION  BY  STAGES  AND  EXAMINATION  OF  RESIDUALS 

(1)  X  REGRINIT  Y 

(2)  REGR  L 

( 3 )  SHOE  V 

(4)  SIRES 

( 5 )  VARIANCE 

' REGRIEIT '  IS  USED  ONCE  AT  THE  OUTSET  TO  SET  UE  GLOBAL  VARIABLES  FOR  'REGR'  . 
'REGR1  PERFORMS  REGRESSION  Oil  ONE  OR  MORE  DESICCATED  INDlISRDEiJT  VARIABLES.  AND 
MAY  BE  CALLED  REPEATEDLY.  'SHOD'  MAY  BE  USED  AT  AtlY  STAGE  TO  OBTAIN  SUMARY  IN¬ 
FORMATION  ABOUT  A  VECTOR.  ' STRES'  GIVES  STANDARDIZED  RESIDUALS  OF  THE  DEPENDENT 
VARIABLE  (S)  FOR  USE  IN  SCATTERPLOTS .  'VARIANCE'  YIELDS  THE  CONVENTIONAL  ESTIinA- 


(1)  RF  GRIN  IT.  THE  FIRST  ARGUMENT  (X)  IS  A  MATRIX  WHOSE  COLUMNS  LIST  VALUES 
OF  THE  INDEPENDENT  ( EXOGENOUS ,  PREDICTOR)  VARIABLES .  USUALLY  ONE  COLUMN  IS  ALL 
VS,  AND  ITS  INDEX  NUMBED  IS  THE  ARGUMENT  IN  THE  FIRST  CALL  OF  'REGR' .  EACH  CO¬ 
LUMN  OF  X  SHOULD  HAVE  BEEN  MULTIPLIED  BY  A  POWER  OF  10  SO  THAT  THE  UNIT  PLACE  IS 
THE  LAST  SIGNIFICANT  ONE.  WHEN  A  COLUMN  OF  X-RESIDUALS  HAS  SUM  OF  SQUARES  LESS 
THAN  (NUm)  IT  IS  JUDGED  TO  BE  ‘TOO  CLOSE  TO  ROUND-OFF  ERROR  TO  BE  WORTH  USING. 

THE  SECOND  ARGUMENT  (Y)  IS  EITHER  A  VECTOR  LISTING  VALUES  OF  ONE  DEPENDENT 
VARIABLE  OR  A  MATRIX  WHOSE  ROWS  LIST  VALUES  OF  SEVERAL  DEPENDENT  VARIABLES  TO  BE 
STUDIED  IN  PARALLEL.  THE  NUMBER  OF  OBSERVATIONS  OR  DATA  SETS  IS  EQUAL  TO 
ltpY  -*-*  itpy 

GLOBAL  VARIABLES  IN  ' REGRINIT '  AND  'REGR': 

P  THE  NUMBER  OF  INDEPENDENT  VARIABLE'S,  (pJO[2]  . 

IND  A  LOGICAL  VECTOR  OF  LENGTH  P  SHOWING  WHICH  INDEPENDENT  VARIABLES  HAVE  BEEN 
BROUGHT  INTO  THE  REGRESSION  RELATION  ( INITIALLY  ALL  0'S). 

HU  THE  NUMBER  OF  RESIDUAL  DEGREES  OF  FREEDOM  ( INITIALLY  THE  NO.  OF  OBS.). 

RX  TRANSFORMED  OR  RESIDUAL  X- VARIABLES  ( INITIALLY  EQUAL  TO  X) . 

RY  RESIDUALS  OF  Y  ( INITIALLY  EQUAL  TO  Y). 

B  REGRESSION  COEFFICIENTS  ( VECTOR  OR  MATRIX)  OF  Y  ON  COLUMNS  OF  X  CONSIDERED 
SO  FAR.  {B  APPEARS  ONLY  IN  THE  OUTPUT  OF  'REGR' .) 

RB  DITTO  WITH  RX  Til  PLACE  OF  X  ( INITIALLY  ALL  0'S). 

TRX  P-BY-P  MATRIX  TRANSFORMING  \  TO  RX.  (qRX  «-*■  THX+.*W  ;  D  RD+.*TRX  . 

THE  FITTED  VALUES  ARE  Y-RY  «-*■  B+.^X  «-»  RB+.*WX  . 

ISS  INITIAL  SUM  OF  SQUARES  (APPEARING  ONLY  IN  THE  OUTPUT  OF  'REGRINIT'). 

RSS  RESIDUAL  SUM  OF  SQUARES  (APPEARING  ONLY  IN  THE  OUTPUT  OF  'REGR'). 

AFTER  'REGRINIT'  IS  EXECUTED  IT  CAN  BE  ERASED,  AS  ALSO  ITS  ARGUMENTS  (WHICH 
SHOULD  BE  STORED  FOR  COPYING  LATER  WHEN  NEEDED). 

(2)  im.GR-  THE  ARGUMENT  (L)  IS  A  SCALAR  OR  VECTOR  LISThC,  THE  INDEX  NO.(S)  OF 
WE  INDEPENDENT  VARIABLE (S)  TO  BE  BROUGHT  NEXT  INTO  THE  REGRESSION. 

RENAME  RSS  OR  ANY  OTHER  GLOBAL  OUTPUT  VARIABLE  TO  SAVE  IT,  BEFORE  AGAIN  CALL¬ 
ING  'REGR'. 

THE  POSSIBLE  REDUCTIONS  Ill  RSS  TABULATED  IN  TilE  1  ‘RUTTED  OUTPUT  SHOW  THE  EF¬ 
FECT  OF  BRINGING  ONE  FURTHER  INDEPENDENT  VARIABLE  INTO  THE  REGRESSION,  A"D  MAY 
BE  USED  TO  GUIDE  THE  NEXT  CALLING  OP  'REGR'. 


THE  SO-CALLED  MODIFIED  GRAM -SCHMIDT  ALGORITHM  IS  USED.  SEE 
AllSCONlBE :  J.  ROYAL  STATIST.  SOC.  B  20  (1967  ),  1-52,  ESPECIALLY  SECTION  1.4. 
BJORCK :  BIT  -  NOR?).  TIDSKR.  INFORMA  T  IONS  DEHAND .  7  (1967),  1-21. 
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IF  YOU  WAtJT  TO  FOIL  TllL  (JU* 4)  PROVISION  XZUTIOi  JU  AdOVi RCPLACL  Tilt'  ,4I  I// 

LINE  Cl]  BY  A  LARGER  NUMBER. 

(3)  SHOW.  THE  ARGUMENT  (l')  IS  A  VECTOR,  SUCH  AS  RY  (IF  A  VECTOR)  OR  A  ROU  OF 

RY  (IF  RY  IS  A  MATRIX)  OR  A  COLUMN  OF  RX  OR  DIAGQ  (SEE  'STRES'). 

NU  IS  USED  AS  DIVISOR  IN  FINDING  THE  R.M.S.  VALUE.  IF  'SHOW'  IS  USED  OUTSIDE 

THIS  REGRESSION  CONTEXT,  HU  MUST  BE  SPECIFIED  FIRST. 

THE  FREQUENCY  DISTRIBUTION  IN  THE  PRINTED  OUTPUT  IS  OVER  6  INTERVALS  OF  EQUAL 

LENGTH.  THE  1ST  AND  6 Tti  INTERVALS  ARL  CENTERED  ON  THE  LEAST  AND  GREATEST  VALUES 

OCCURRING  IN  V. 

(4)  STRES.  NO  ARGUMENTS ;  SHOULD  BE  USED  ONLY  AFTER  'REGR'  HAS  BEEN  EXECUTED. 

GLOBAL  OUTPUT  VARIABLES : 

DIAGQ  IS  A  VECTOR  CONSISTING  OF  THE  MAIN  DIAGONAL  OF  A  MTRIX  Q  (NOT  FURTHER  DE¬ 
TERMINED  HERE)  THAT  PROJECTS  Y  INTO  RY,  THUS  RY  Y+.*Q  . 

SRY  EACH  ELEMENT  OF  RY  IS  DIVIDED  BY  THE  SQUARE  ROOT  OF  THE  CORRESPONDING  ELE¬ 
MENT  OF  DIAGQ  (TIMES  A  CONSTANT ).  IF  THE  MEMBERS  OF  DIAGQ  ARE  ALL  EQUAL, 

SRY  IS  THE  SAME  AS  RY. 

(5)  VARIANCE.  NO  ARGUMENTS ;  ISE  ONLY  AFTER  'REGR'  HAS  BEEN  EXECUTED. 

ONE  GLOBAL  OUTPUT  VARIABLE: 

XPXI  IS  ' X-PFIME  X  INVERSE' .  UHEN  REGRESSION  HAS  BEEN  PERFORMED  ON  ALL  COLUMNS 
OF  X,  XPXI  IS  THE  INVERSE  OF  ($ *)+.**;  OTHERWISE  IT  IS  THE  CORRESPON¬ 
DING  EXPRESSION  FOR  THE  COLUMNS  OF  X  THAT  HAVE  BEEN  USED,  TOGETHER  l/ITH 
ZEROS  FOR  REGRESSION  COEFFICIENTS  NOT  ESTIMATED. 
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V  X  REGRIN IT  Y\J 

[1]  ->3-(A/iSpX)A(A/ispy)A(v/  l  2  =ppY)*2=ppX 

[2]  ->3+J+-(~l*pY)=NU+lipX 

[3]  -^0,p|>'/70  GO. ' 

[  4  ]  X+Oi,  TRX+(P  ,P  )p  1  ,INEh-(  P*-~  1  ipRX+X )  p  0 

[5]  -*-9xxi-pPiss<-+//?yx/?y*-y 

[6]  RB*-PpO 

[7]  Q,' INITIAL  SUM  OF  SQUARES  (ISS),  NUMBER  OF  READINGS  AND  MEAN  SQUARE  ARE' 

[8]  -*-0,i>(L0.5+IS5),/W,0.001xC0.5+1000xJ,S5T//y 

[9]  RBH(UpY),P)pO 

[10]  Q,'  *Y-VARIADLC  J 

[11]  'INITIAL  SUM  OF  SQUARES  ( ISS(  NUMBER  OF  READINGS  AND  MEAN  SQUARE  A 

RE' 

C 1 2 ]  ( l0. 5+JS5CJ]  )  ,/W ,0.001*  LO .  5  +  1000 xJ55[(/]t/W 

[13]  -*-iox(itpy)^j+i 
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7  STRES ’,M 

[1  ]  ->2+ a/0  <DIAGQ*-\—t-  //•/:-  ( pM)  p+//-M  IND/RX  )*2 

[2]  -*-0,pOM  MEMBER  OF  DIAGQ  VANISHES .  SRY  NOT  FOUND.  ' 

[  3  ]  SRY+-RY*  (pRY)p(DIAGQ*(pRX)ll]HJU)*~0.5 
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?  REGR  LiDiliJ-JiiSiU 
Cl]  -»2+(0<i><-4)Aj-*-(A//Je\P)AA/lSpiA,Z, 

[2]  •♦O.pLh-W  GO.  ' 

[3]  -4+1  <,NU+-NU- 0  =IND[  K+-LIJ  ]  1 
[  4  ]  -*-0 ,  pi> '  NU  HAS  VANISHED .  ' 

[5]  awin'}*  1 

[6]  -*9*  i  (//[/*D)sS^r  ;\']t .  *RXi  \K] 

[7]  • DIVISOR  TOO  SMALL  AT  COL  Urn  1  '  ;K;'  OF  RX. ' 

[8]  -0 

[9]  RX+RX-RXl  ;A']°  ;A']k  *RX*{pRX)f}~IIID)*5 

[10]  raMHH/»,»Tffl[i:;] 

[11]  -12+2=pp/?Y-/?Y-((M/?Y+.x/?A'[  ;A'])*S)°.xtfY[;A] 

[12]  -14,/?PU>®[A']+U 

[13]  RBlua+RDliKhU 

[14]  +3*\(pL)ZJ*-J+l 

[15]  I<-(NUiDHS++SHX*RX 

[16]  RSS<-+/RY*RY 

[17]  -22x \J*-v/2=ppd+fio  + .*TRX 

[18]  Q,' RESIDUAL  SUM  OF  SQUARES  (RES),  D.F.  ( tJU )  AND  NEAR  SQUARE  ARE ' 

[19]  ( [0. 5+ASS)  ,W/,O.OOlx[o ,  5+1000*/?SS ri.U 

[20]  • POSSIBLE  REDUCTIONS  IN  RSS  ARE ' 

[21]  -28 ,1>L0. 5+((Jx/?Yt .  *RX)*  2 )  *S[NU*D 

[22]  *Y •VARIABLE  ' 

[23]  ' RESIDUAL  SUM  OF  SQUARES  (RSS[_ 0.F.  (M7)  AND  MEAN  SQUARE  ARE' 

[  24  ]  ( L  0 . 5+/?5S[«7  ] ) ,  HU ,  0 . 0  0 1  x  L  0 . 5+ 1 0  00  xRSSlJ]  U1U 

[25]  'POSSIBLE  REDUCTIONS  Ill  RSS{  '  ; J ; '  ]  ARE ' 

[26]  IQ. 5+({I*RY[J  ib  .*RX)*2)iS{NUiD 

[27]  -*-22xi(l+p/?y)>(7^f/+i 

[28]  -Oxil=A  n 

[29]  Q, ' DIVISOR  TOO  SMALL  AT  COLUMN  '  ;(~I)/iP;'  OF  RX.  ' 
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[1] 

[2] 

[3] 

[4] 

[5] 

[6] 
[7] 

[B] 


V  SHOW  V; S\ R;l/ 

■*3“  (v/6SpI7)Alrppl/ 

-3+0<-M(/?-UlO[(i3),(l+plO-<tn3])[6  1]] 

-0,pO 'NO  GO.  ' 

S-K).01xL0.5+100x((+/yxlO*M/)*0.5 

^•K).01xL0.5+100xH;/f] 

•3  LOWEST,  3  HIGHEST  {WITH  RANKS)  AND  R.M.S.  VALUES' 

Iff/;’  ('iltfli'),  '  ;i/[2]  ; '  (  '  ;/f[  2  ] ; ! ) ,  ’;l/[3];'  C;/f[3];'),  '  ;f/[4]  ;  *  (';/?[ 
4];').  ';W[5];'  (';/?[  5];')  ( ’ ;  5  +/?;•);  '  ;5 

'FREQUENCY  DISTRIBUTION :  '  ; +  /( 0 , 1 5 )°  .  =1 0.5+D*(^-V[ltff])*-/W:«[6 

1]] 

V 


V  VARIANCE 

[1]  YPXI-(  ( (P,P)oIND\i+/IND/RXxRX)*<SlTRX)  +  .  *TRX 

[2]  -(12  =pp/?Y)/  3  4 

[3]  -+Q,pLh-' ESTIMATED  VARIANCE  MATRIX  OF  J  IS:  XEXI*RSSMW' 

[4]  'ESTIMATED  VARIANCE  MATRIX  OF  DUih  IS:  XPXI*RSSlJ]*NU' 

V 
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V  RONCOL  y  ;N;A 

Ll  J  ‘♦2+(bS*/py  )A(A/2s|)i’)Ai;  =  (.py 

1 2 ]  C().' 

[3]  £,  'GRAND  MLA.<  (GM)  IS  '  ;GM*-{  +  /,YK.  ~*/pl' 

[4]  ' ?£>m  Ji/w  OF  SQUARES  AvOUT  MEAN  (TSS),  DEGREES  OF  FREEDOM  Aid)  ML'Au  SQUARE 

ARE  '  -,TSS,(N-\ ) ,(  i/V-1  )*TSS-+/ ,RY*RY*-Y-Gt1 
C  53  Q, '  */iW5*’ 

[  b  ]  '  EFFECTS  ( /?£' )  A/fz:  •  ; /?£'-(  +  /RY )  *  H / 

[7]  'SUM  OF  SQUARES  (SSR),  U.F.  hiJDM.S.  ARE  '  ;55/?,  Cl-UlpY ) ,  (  ri  +  ltpy)*55fl«-(  1 
+pY  )*/?£'+ .  »/?£ 

Ld]  Q, '  *  COLUMNS* ' 

[<J]  'EFFECTS  (CE)  ARE  '  iCtX+tRY  )iUpY 

[10]  'SUN  OF  SQUARES  (550,  D.F.  AND  M.S.  ARE  '  ;SSC ,(~  1+UpY)  ,(i~UUpY)*SSC*-(l 
tp Y)*CE+.*CE 

[11]  £,’  * RESIDUALS*' 

[12]  'SUM  OF  SQUARES  (RSS),  D.F.  (NU)  AND  M.S.  ARE  '  -,RS5,t’U,  ( *FUW”l+py)*fl5S«-+ 
/  ,RY*RY*-RY-RE° .  +05 

[13]  '77/F  Ml  77?/ A  OF  RESIDUALS  IS  SAVED  RY.  ' 

[14]  -K0<55tfx55C)/_2+pLK:, '  *TUFJ'/"S  TEST*' 

[15]  +&}*■' TEST  FAILS  BECAUSE  SSROSC  VAN  I  SUES. ' 

[lb]  'REGRESS  I OU  CQF.FFT.  (o')  OF  57  0/.'  RF*.*CR  IS  '  ;/M /!«-+/  ,RY*RF°  .  *CL)*tMSSR* 
SSC 

[17]  'TV KEY'  '5  1-D.F .  TERM  IS  ' 

[lb]  'REMAINING  SUt-  OF  SQUARES,  D.F.  AND  M.S.  ARE  '  ;/! .  (,70-1 ) ,  (  :  NU-1)*A+-RSS-A 
[19]  a  iTf//5  PROGRAM  PERFORMS  A  STANDARD  A.  A  LYSIS  OF  VARIANCE  ON  A  TOO-' 'AY  TAEL 
E. 

V 


V  ROUCOLPERMUTE 

[1]  RE+RSlRMRb'} 

[2]  CE+CElCP+lCE] 

[3]  RY+RYIRP-.CP] 

[4]  'ROW  PERMUTATION  ( RP )  IS  '  ;RP 

[5]  'COLUMN  PERMUTATION  (CP)  IS  ' ;C T 
[b]  'RE,  CE  AND  RY  ARE  NON  PERMUTED .  ’ 

[7]  A  USE  AFTER  ROM  COL,  DEFURL  ROMCOLDI  SPLAY . 

V 


V  ROWCOLDI SPLAY  i\C\D-,R\S\i: 

[1]  >2+(0<t/,J)A(0;ppi)A(2=p5-  1  l)A(/)*-2ki2 

[2]  -*>0,pO 'NO  GO.' 

[3]  -*4  +  13O2ptV("itf/«-(D,0*ir[y.‘j+((  UO;')-"HO'.’)t/)+  .»(  ipOZ.>.  <ipCDpO 

[4]  -*-p|> 'THE  ARGUMENT  IS  TOO  SMALL,  TRY  AGAIi  .' 

[5]  C[0>  1 

[b]  *CODED  DISPLAY  OF  PERMUT'D  RES  I  DU  ALU*' 

[7]  RH  ltlMl,ir[0.r.+((‘HFF)-H/?/’)n//,F)t.x(  xpRL')«  .  <i  pFF)pO 

[8]  F[W]*-1 

[y]  'NE-SU  DISPLACEMENTS  ROUGNLY  nL'ASURE  CHANGES  IP  FITTED  VALUES. '  ,Q 

[10]  CV?\,tf/-°+P£,[[4.St3*Fy;r/,|/i7] 

[11]  A  SEE  '  tlOUTRIPLEPLOT ' . 

V 


1 


V  ANALYZE  Y 

[1]  -*-2+(A/2spy)A2<ppy 

[2]  -*0,pO 'NO  GO .' 

[3]  w«-(i,itppy)pi 

[4]  'GRAND  MAN  ( CM )  IS  '  ;aM  +  /,y)  WpY 

[5]  TOm  SUM  OF  SQUARES  ABOUT  MEAN  ( TSS ),  DEGREES  OF  FREEDOM  AND  MEAD  SQUARE 

ARE  '  iTSS,UUA*M*-~U*/pY)*TSS4-+/,RY*RY-Y-GM 

[6]  'PROCEED  DY  REPEATEDLY  CALLING  THE  FUNCTION  "EFFECT".' 

[7]  'THE  ARRAY  OF  RESIDUALS  IS  ALi/AYS  NAMED  RY .  1'ITH  DEGREES  OF  FREEDOM  NU 

[8]  'RENAME  SS  EACH  TIME  TO  SAVE  IT.  ' 

[y]  n  ZIEGI//S  4a’  ANALYSIS  OF  VARIANCE  OF  A  PERFECT  RECTANGULAR  ARRAY. 

V 


V  Z+EFFECT  V \J \K\ii\P\R 

[1]  ->3-a/v/M*-(  ,V)°.  =  \R++/ppRY 

f2]  -*-3tA/0</l>*-l-+jW 

[3]  -O.p  Q-'NO  GO.' t2+" 

[4]  J*-t/ltQDF+DF,[l)  M, 0 

Cb]  SS+-R**/ ,Z*ZHiK)x*/((K+*/M/pRY),(pRY)[V])p(W<-(M/\R),V)X)RY 

c  7  ]  RY+RY-m  ( w/p/?y ) ,  (P/?y )  c  vi  >Pz 

[8]  'SUM  OF  SQUARES  (SS),  D.F.  ANDM.S.  ARL  ' -,SS,DFU ;UR),SSiDFlJ;l+R} 

[9]  a  • ANALYZE '  SHOULD  BE  CALLED  FIRST. 

V 


V  Y*STDIZE  X-.C 
Cl] 

[2]  -K),pl>'/70  GO.'.y-*-" 

[3]  y+y*((T~i+Px)x+/yxy+x-(+/x)*pA)*''o.5 

[4]  a  27/£  VECTOR  X  IS  RESCALED  TO  HAVE  ZERO  MEAN  AND  UNIT  VARIANCE. 

V 


\ 

\ 


HOWFILTER 


FILTERING  TIME  SERIES 

(1)  IM1  FILTER  X 

(2)  (Ml  MV  X 

(1)  FILTER .  TtlE  FIRST  ARGUMENT  (,/)  IS  A  PAIR  OF  INTEGERS  DEFINING  THE  FILTER. 
THE  SECOND  ARGUMENT  ( X )  IS  A  VECTOR  OF  DATA  TO  RE  FILTERED.  FILTERING  CONSISTS 
OF  SUBTRACTING  A  COSINE-WEIGHTED  MOVING  AVERAGE  OF  EXTENT  Ml 2]  FRO,,  A  SIMILAR 
MOVING  AVERAGE  OF  EXTENT  A/C  1  ]  .  THE  RESULTING  HEIGHTS  ARE  DISPLAYED,  TOGETHER 
WITH  THEIR  SUM  uF  SQUARES.  THE  ELEMENTS  OF  U  MUST  HE  EITHER  DOTE  ODD  OR  BOTH 
EVEN,  OR  ELSE  ONE  OF  THEM  MUST  UE  0.  TO  GET  Ail  11  -POINT  MOVING  AVERAGE.  SET  M 
EQUAL  TO  11  0  .  TO  DETREND  A  SERIES  BY  SUBTRACTING  FROM  IT  AN  11-POINT  MOVING 
AVERAGE ,  LET  M  BE  1  11  .  FOR  A  SMOOTHED  DETRENDED  SERIES ,  LET  M  BE  (E.G.) 
3  11  .  SETTING  M  EQUAL  TO  3  1  YIELDS  THE  SECOND  DIFFERENCE  OF  X  ( DIVIDED 
BY  4).  IF  THE  HON  ZERO  ELEMENTS  OF  M  ARE  ODD,  THE  RESULT  CORRESPONDS  TO  THE 
SAME  TIME  POINTS  AS  THE  DATA  ( APART  FROM  LOSS  OF  SOME  END  VALUES).  BUT  IF  DOTH 
ELEMENTS  OF  M  ARE  EVEN ,  THE  RESULT  CORRESPONDS  TO  TIME  POINTS  MID, 'AY  BETWEEN 
THOSE  OF  THE  DATA. 

IF  AN  ELEMENT  OF  M  EXCEEDS  4 ,  A  LITTLE  CHEATING  GOES  ON  AT  THE  ENDS  OF  X  , 
THE  FIRST  AND  LAST  ELEMENTS  BEING  REPEATED  R  TIMES  AS  DEFINED  AT  [7].  IF  YOU 
DO  NOT  NANT  THIS  FEATURE,  EITHER  REPLACE  [7]  BY 
[7]  R-*-Q 

OR  DROP  THE  FIRST  AND  I  AST  R  ELEMENTS  OF  U  . 

(2)  Aj&K.  AN  ADAPTATION  OF  ' FILTER •  TO  TAKING  A  COSINE-NEIGUTED  MOVING  AVERAGE 
OF  EACH  OF  SEVERAL  VECTORS  SIMULTANEOUSLY .  THE  FIRST  ARGUMENT  (M)  IS  A  POSITIVE 
INTEGER,  THE  SECOND  ( X )  IS  A  MATRIX.  THE  RESULT  (U)  IS  A  MATRIX,  WHOSE  RONS  ARE 
Mi  AT  WOULD  BE  OBTAINED  BY  REPEATED  APPLICATION  OF  ' FILTER *  ( MODIFIED  AS  SUGGEST¬ 
ED  ABOVE  BY  SETTING  R+ 0)  TO  THE  ROWS  OF  X  : 

UlJ’,1  ~  (A/,0)  FILTER  X[J ;  3 

IF  YOU  WANT  THE  HEIGHTS  TO  BE  DISPLAYED,  INSERT  +  AT  THE  BEGINNING  OF  LINE  [4] 
OR  SOMETHING. 

28  JUNE  1973 


V  (Ml  FILTER  X;P,R;W 

[1]  -*-3-(A/2=p/'/)A(i=ppA/)Al=ppX 

[2]  +3W(0SM)A(A/=lA/)A((0  =  K/A/)v0=2|-//‘/)A(*/A/)A(2*~2+pJO2P<-r  /M 

[3]  -Hl.pO  ,NOGO.',U+" 

[4]  /ML0.5xP-A{l])<K(P-A/tl])p0),<l-2o(iWCl])xO2U+A/tl])rl+W[l] 

[  0  3  l/Ht- (10. 5  *P-Af[  2  ]  )4>(  ( P-Ml  2  ]  ) p 0  )  ,  (  1  - 20(  t  A/[  2  3 ) xo2  r  1 +M[  2  J  H 1  +AJ[  2  3 

[03  'FILTER  WEIGHTS  ARE  '  ;i/;  '  WITH  S.S.  x,7 

[73  /?H(F-1H4 

[83  ■+£tl*i(l22)<1b*P*p/Mflpitx)  ,X,Rp~lfX 
[9]  +0,p(M/+.x(0,l-P)4.Cl+iP)(|>(P,pA'>pX 

[10 J  £l:fl«-pp(M(pA)+l-P)pO 
[113  (A*W[fl3*(tf-i)+(ff-P)+X 
[12]  HP2R+R+l)/Ll+l 

7 
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UOWHARMONIC 

FOURIER  ANALYSIS  OF  A  TIME  SERIES 

(1)  V+K  TAPER  U 

( 2 )  Z+ilA  V 

(1)  TAPER.  THE  FIRST  ARGUMENT  (K)  IS  A  POSITIVE  INTEGER,  AND  THE  SECOND  (U) 
IS  A  VECTOR  OF  LENGTH  GREATER  THAN  2*K  .  'HIE  RESULT  O')  IS  A  VECTOR  OF  LENGTH 
(pU)-K  ,  A  CIRCULARIZED  VERSION  OF  U  OBTh+i  ED  NY  OVERLAPPING  THE  FIRST  AND  THE 
LAST  K  ELEMENTS.  THE  FIRST  K  ELEMENTS  OF  U  MULTIPLIED  BY  A  VECTOR  OF  LI¬ 
NEAR  HEIGHTS  N  ,  PLUS  THE  LAST  K  ELEMENTS  OF  U  MULTIPLIED  BY  1-W  ,  FORI!  K 
ELEMENTS  OF  V  .  SOME  AT  THE  BEGINNING  AND  SOME  AT  THE  END.  THE  REST  OF  V  IS 
THE  SAME  AS  THE  REST  OF  U  .  THE  WEIGHTS  V  ARE  DISPLAYED.  DIFFERENT  WEIGHTS 
CAN  BE  OBTAINED  BY  ALTERING  THE  DEFINITION  IN  LINE  [4].  THIS  PROCESS  IS  NEARLY 
BUT  NOT  QUITE  WHAT  IS  GENERALLY  CALLED  TAPERING.  TO  PERFORM  A  HARMONIC  ANALYSIS 
OF  A  GIVEN  TIME  SERIES,  ONE  MAY  SUCCESSIVELY  EXECUTE  'FILTER' ,  'TAPER'  AND  'HA'. 

(2)  dA-  THE  ARGUMENT  ( V )  IS  A  VECTOR  OP  LENGTH  N  NOT  LESS  THAN  3.  THE  RE¬ 
SULT  (Z)  IS  A  MATRIX  WITH  TWO  ROWS.  THE  FIRST  ROW  IS  THE  FOURIER  DECOMPOSITION 
OF  THE  TOTAL  SUM  OF  SQUARES  OF  V  ( ESSENTIALLY  THE  LINE  SPECTRUM  OR  PERIODOGRAM) . 
THE  SECOND  ROW  LISTS  PHASES.  THE  ANALYSIS  OF  VARIANCE  SHOULD  BE  CHECKED  FOR  NU¬ 
MERICAL  ACCURACY: 

+/Z[l ;]  «-*  V+.*V 

ZClil]  IS  THE  'CORRECTION  FOR  THE  MEAN •  WITH  1  DEGREE  OF  FREEDOM.  FOR  NUME¬ 
RICAL  STABILITY  IT  OUGHT  TO  BE  SMALL.  THE  CORRESPONDING  PHASE  Z[2;l]  IS  SET 
CONVENTIONALLY  EQUAL  TO  0.  IF  N  IS  EVEN,  THE  LAST  MEMBER  OF  Z[l;]  ALSO  HAS 
1  D.F.  AND  THE  CORRESPONDING  PHASE  “l^Zt2;3  IS  SET  CONVENTIONALLY  EQUAL  TO  0. 
ALL  OTHER  MEMBERS  OF  Z[l;]  HAVE  2  D.F.  Z[l;J+l]  REFERS  TO  THE  J  -TH  HARMO¬ 
NIC,  J  CYCLES  PER  N  TIME  UNITS,  WHERE  J  IS  LESS  THAN  Hi 2  .  PHASE  ( BETW¬ 

EEN  0  AND  02)  IS  MEASURED  RELATIVE  TO  A  TIME  ORIGIN  ON  THE  LEFT,  SO  THAT  THE  OB¬ 
SERVED  TIME  POINTS  APE  u)  .  THE  FITTED  HARMONIC  IS 
AMPLITUDE  x  COSINE  ( PHASE  +  (02*JiN)*\N) 

WHERE 

AMPLITUDE  *-  ( Z[  1  ;</+ 1  ] x  2  f N ) * 0 . 5 
PHASE  «-  Z[2 ;J+l] 

15  NO V.  1971 


Z+HA  V;J \M\N\A 

-*2+(  3  $//•*-+ /p  10  a  («/«-l  )=ppl/ 

-*0,pL B-'NO  GO.'  ,Z+" 
ZH2,UUP'2)p((A*A++/V)lN),(N+l)pO 
Fh-  1  2  °.0(02iN)*u! 

ZCl  ;J+lX2*//)x/k.x/W;[  ;i+//|”i+Jx u;]+.xy 
-*8-0=/l[2] 

-•■9 ,  Z  [  2 ;«/+ 1  >oo .  5+ 1 = x/1  [  1  ] 

Z[2;</+l]«-(02)  j  (o"l=x/i[2]  )-~30*//l 
-*(10  1  =*/i/-2xJV+l)/  5  10  0 

ZCl  ;t/+l  ]+{A*A*-/V)iN 


10  - 


- - r-v  - - ■  -■  ,  ■-  TJfla--.: 


V  IM4  MAV  X-,R\W 
Cl]  +3-(0=ppM)A2=ppX 

[2]  -*>3+(W21)A(/./=LW)AAi£“ltpX 

[3]  -KJ.pLH-'/lAF  NO  GO. '  ,(/«•'  • 

[4]  O(l-2o(  i;/)xo2*l+A/)fl+M 

[5]  Mi '  -POINT  COSINE -WEIGHTED  MOVING  AVERAGES1 

[6]  OpppO(  (pA')-OfM-l  )pO 

C7]  £l:OOfc/[ff]x(o,/?-l)i(0,/?-M)+X 
[8]  -+(a/2/W?+1)/£1 

V 


V  OK  TAPER  U;(/ 

[1]  -KM0=pp;OAl=ppi/ 

[2]  -*-3+(A/=U{)A(^l)Ali(pi/)-2xX 

[3]  -*0,p[ Tr'tiO  GO.' tV+" 

[ 4 ]  '  TAPER  WEIGHTS  ARE  ';(Mi A' )*X+1 

[5]  //[iK>((-K)tA)+f/x//[iK]-(-A')+/; 

C6]  0<W*2  )$(-*)«/ 

V 


V  0/1/  AUTOCOV  V-,J 

[1]  ‘►2+0/£"l++/pVr)A(O2£+/pl0A(l-pp’/)A0=pp/i/ 

[2]  -K),pl>'//0  GO.  ’.O" 

[3]  or+.xy 

[4]  -*4x/teOpOC,(Jil')+.>‘(-</)  +  l' 

[5]  p  FIRST  H  TERMS  OF  CORRL'LOGRAM  OF  V  ARE  HOC[l] 

V 


V  Z*-//!?  X 

[1]  Z«-0.5+(*X)x0.5-(*X*Xr2)x((U+0.23164l9*|X)°.M5)+.x  319381530  ”356503782 
1781477937  ”1821255978  1330274429  x3 . 989422804S”l0 

[2]  n  HASTINGS  APPROXIMATION  TO  THE  NORMAL  {GAUSS -LAPLACE)  INTEGRAL  FUNCTION , 

AftS  HANDBOOK ,  /WS55,  26.2.17. 

[3]  P  ABSOLUTE  ERROR  IN  Z  LASS  7YMA  1£~7.  X  AM 7  BE  ANY  NUMERICAL  ARRAY. 

V 


V  X+-INIF  PiU 

[1]  ->2+  a/  ,  ( 1  >P )  ao  <P 

[2]  +pl>,«0  G.9. ' 

[3]  O(~2*»0. 5  -  |P  *0 5  )*0. 5 

[4]  X-*-(xP-0.5)x//-((A°.*'l+i3)+.x  2515517  802853  1032b )i{U»  .*~l+i4  )+.* 

1000000  1432788  189260  1308 

[5]  P  HASTINGS  APPROXIMATION  TO  THE  INVERSE  OF  THE  NORMAL  INTEGRAL  FUNCTION , 
NBS  HANDBOOK,  AMS 55,  26.2.23. 

C6]  p  ABSOLUTE  ERROR  IN  X  LESS  THAN  5A~4.  P  MAY  BE  ANY  ARRAY  OF  NUMBERS  BETW 
EEN  0  AND  1. 

V 


ll 


i, mi  «,..l  I, up  l.ii.  i 


^5»ww,  wwii,i7  ■«  j^*i  '  nr- 


tiOWHUBER 

ROBUST  REGRESSION 

(1)  Z Z+R  HUBER  Z 

( 2 )  DZ+HUBERl 

( 3 )  DZ+-HUBCR2 

THE  FUNCTION  'tlUBER'  PERFORMS  ORE  CYCLE  OF  ITERATION’  TOUARDS  MINIMIZING  THE 
SUM  OF  A  FUNCTION  RHO  OF  THE  RESIDUAL L  IE  A  REGRESS IOE  PROBLEM ,  l /HERE  RHO  IS 
DEFINED  IN  TERMS  OF  A  POSITIVE  CONSTANT  K  AS  FOLLOl/S : 

V  U  <-  RtlO  Z 

Cl]  fM(0.S*Z*2W<2|Z)  +  (A,><( \Z)-Ki2)*K<\Z  V 

(SEEP.  J.  HUBER:  &UN.  AM  77/.  STATIST.  35  (1964),  73-101.) 

K  (MY  BE  SPECIFIED  IN  ADVANCE,  OR  ALTERNATIVELY  MAY  BE  DETERMINED  BY  THE  CON¬ 
DITION  THAT  A  SPECIFIED  PROPORTION  OF  RESIDUALS  EXCEED  K  IN  MAGNITUDE.  CALCU¬ 
LATION  OF  ESTIMATED  CHANGES  IN  THE  REGRESSION  PARAMETERS  AND  IN  THE  RESIDUALS  IS 
DONE  BY  A  SUBSIDIARY  FUNCTION ,  'HUBER!'  OR  'HUBERT1.  IF  NEITHER  OF  THESE  FUNC¬ 
TIONS  IS  APPROPRIATE  FOR  THE  DATA,  INSERT  A  LIEU  LINE  IN  ' HUBER '  AS  FOILS* IS : 

[8.1]  +44.P DZ+HUDER2 

AND  DEFINE  A  NEW  SUBSIDIARY  FUNCTION  STARTING 

V  DZHJUBER2 ;  ...  V 

(1)  HUBER.  THE  Fa.  IT  ARGUMENT  (R)  MUST  BE  SCALAR.  IF  IT  IS  ~1,  THE  CURRENT 
GLOBAL  VALUE  FOR  K  IS  USED-,  OTHERWISE  R  MUST  LIE  BETWEEN  0  AND  1,  AND  IS  THE 
PROPORTION  OF  RESIDUALS  TO  EXCEED  K  IN  MAGNITUDE.  THE  SECOND  ARGUMENT  (Z)  IS 
THE  ARRAY  OF  RESIDUALS  CORRESPONDING  TO  A  TRIAL  SETTING  OF  THE  REGRESSION  PARA¬ 
METERS  BETA  ,  WHICH  MUST  BE  SPECIFIED  BEFORE  ' HUBER'  IS  CALLED.  THE  EXPLICIT 
RESULT  (ZZ)  IS  A  NEW  SET  OF  RESIDUALS,  THAT  CAN  BE  USED  AS  SECOND  ARGUMENT  IN 
THE  NEXT  CALL  OF  'HUBER' . 

GLOBAL  OUTPUT  VARIABLES:  K,  DUET  A  (ESTIMATED  REQUIRED  CHANGE  IN  PARAMETERS), 
BETA  (NEW  PARAMETER  VALUES). 

THE  PARAMETER  CHANGE  DBETA  YIELDED  BY  THE  SUBSIDIARY  FUNCTION  IS  DETERMINED 
FROM  THE  FIRST  AND  SECOND  DERIVATIVES  OF  THE  SUM  OF  RHOS  AT  THE  INITIAL  BETA  . 
IT  YIELDS  THE  TRUE  MINIMUM  OF  THE  SUM  OF  RHOS,  FOR  THE  VALUE  OF  K  USED,  IF  NO 
CHANGE  OCCURS  IN  THE  INDEXING  OF  THE  RESIDUALS  THAT  EXCEED  K  IN  MAGNITUDE  — 
THOSE  RESIDUALS  BEING  REFERRED  TO  AS  'MODIFIED  RESIDUALS' .  OTHERWISE  THE  TRUE 
MINIMUM  IS  NOT  OBTAINED.  THEN  THE  RESIDUALS  ARE  CALCULATED  WITH  THE  FULL  PARA¬ 
METER  CHANGE  DBETA  AND  ALSO  WITH  CHANGES  THAT  ARE  ONLY  0.9,  0.8,...  OF  DBETA, 
STOPPING  WHEN  THE  SUM  OF  RHOS  IS  LEAST.  THE  MULTIPLIER  OF  DBETA  IS  NAMED  AL¬ 
PHA  IN  THE  DISPLAYED  OUTPUT.  THE  VALUE  ALPHA  =0.0  CORRESPONDS  TO  THE  INPUT 
PARAMETER  SETTING  BETA  .  WHAT  IS  FINALLY  YIELDED  AS  THE  OUTPUT  VARIABLE  DBETA 
IS  THE  FIRST  DBETA  MULTIPLIED  BY  THE  BEST  ALPHA  ;  AND  THAT  DBETA  IS  ADDED 
TO  THE  INPUT  BETA  TO  MAKE  THE  OUTPUT  BETA  . 

'  HUBER  V  AND  '  HUBER2 '  WILL  FAIL  IN  EXECUTION  IF  THE  RIGHT  ARGUMENT  OF  ffl  IS 
SINGULAR.  IF  THAI  HAPPENS,  (CLEAR  THE  STATE  INDICATOR  AND)  START  'HUBER'  AGAIN 
WITH  A  LARGER  K  (OR  SMALLER  POSITIVE  R  ). 

(2)  HUBER!  IS  INVOKED  WHEN  Z  IS  A  VECTOR.  THERE  MUST  BE  A  GLOBAL  MATRIX 
X  WHOSE  COLUMNS  ARE  THE  INDEiENVbin  VARIABLES  IN  A  REGRlSsTOn.  TOR  CONSISTENCY 

pZ  ♦->  lip*  ;  p BETA  ++  ~ltpJf  . 

THE  ORIGINAL  DATA  VECTOR,  NOT  NAMED  IN  THE  PROGRAM,  WAS  EQUAL  INITIALLY  TO 
Z  +  X+.xDETA  . 


TiMi 


(3)  HUBER 2  IS  INVOKED  UUCu  Z  IS  A  MATRIX.  2  MUST  HAVE  AT  LEAST  3  flC/S  AI'D 
3  COLUMNS .  THE  USUAL  ADDITIVE  STRUCTURE  IS  FITTED.  BETA  IS  A  VECTOR  OF  LEE  GTE 
l++/p2  ,  AND  LISTS  THE  MEAN,  THE  ROD  CONSTANTS  ( SUmiNG  TO  0),  4/Y?  7V/Z7  COMM/' 
CONSTANTS  (SUMMING  TO  0).  THE  ORIGINAL  DATA  MATRIX,  NOT  NAMED  IN  WE  PROGRAM, 
NAS  EQUAL  INITIALLY  TO 

l  +  Om[l]  +  ((pZ)[l]tHPm)<>.  +  (pZ)[l]mom  . 

23  JiM,’  1973 


V  ZZ+fl  HUBER  Z\DZ\I',I1  ;*/ ;<¥;P;5;i>l;ZZl 

Cl]  -+2+  ( ppOfiT/H-  .OaT/I )  a  ( 0 =p  pF )  AJ+-3  ZN+x/pZ 

[  2  ]  *&1 ,  pO '  i/£/F£ff  NO  GO . ' 

[3]  ->(P='l)/6 

[4]  ■+(  ( l>P)v/,/<p«-0. 5+A*l-F)/2 

[5]  K+~+/  ( 1  Ot‘11  *P-LP)*|(,Z)[UI  ,Z)[(LP)JP]] 

[6]  'A'  =  ';F 

[7]  'INITIAL  NUMBER  OF  MODIFIED  RESIDUALS  =  ';+/,I«-F< |Z 

C8]  'FOR  ALPHA  =  0.0,  SUM  OF  RHOS  =  '  ;(0. 5x+/,Z>«Zx~J)+/fx+/,Ix(  |Z)-Ft2 
C9]  -►(  (1  2  =ppZ)/£2,&3),2 

[10]  4l:-K) ,ZZ«-" 

[11]  42 :  -*>(  41  x  i  o = x  / p DZ-HUBER 1 ) , £4 

[12]  lL3:+{ilx\0=x/pDZ+-HUNER2 

[13]  £41-^46X1  a/  ,J=J1^A'<  |  ZZ+Z+DZ 

[14]  'FOP  ALPHA  =  1.0,  SO//  OF  RHOS  =  '  ;S«-(0 . 5x+/ ,ZZxZ2x~Il  )+A'x+/ , Jlx(  |  ZZ)-Ft 
2 

[15]  45:Il^<|ZZl^Z+OZx^«7-0.1 

[16]  'FOR  ALPHA  =  SUM  OF  RHOS  =  •  ;51*(0.5x+/,ZZl*ZZlxrn)+Kx+/,Il*(  |ZZ1) 

-.{*2 

[17]  ^7xi(OsSl)vJ=0 

[18]  S+-S1 

[19]  -»45,pZZ^ZZl 

[20]  46 :  'MINIMIZED  SUM  OF  RHOS  =  '  ;(0.5x+/  ,ZZxZZx~J)+Fx+/,Jx(  |ZZ)-K*2 

[21]  47 :'DDETA  =  ' ; DDETA+DBETA * 1 L J+0 . 1 

[22]  'Pm  =  '  lUETA+BETA+DiShTA 

V 


V  DZ+HUDER1  ;Y 

[1]  ■♦2+((pPm)="ltp/)A(P=ltpA)A(2=ppZ)Al=ppZ 

[2]  -►pO'/M.'PPl  AO  G0.',DZ<-" 

[  3  ]  0ZW+ .  xDBCTAH  ( A'l  Zr  -A )  + .  xX  )$(  )  + .  x  Y+  ( ~I )  +X 

V 


V  DZ+HUBER2  -,A-,D-,C\D\M-,U 

[1]  *2+((pO£T/l)=l-M7pZ)A(A/3SpZ)A2=ppZ 

[2]  +pQ<-'HUUER2  NO  GO. '  ,PZ<-' ' 

[3]  iM-(+/,y),(+/i/),+/0«-A'lZr-F 

[4]  CM(pZ)[2  2]pl,(pZ)[2]p0)x(pZ)[2  2]p/l^~I 

[5]  ZM( pZ)[l  l]pl,(pZ)[l]pO)x(pZ)[l  l]pP*M-/~J 

[6]  /*M(+/P). -(+/!),+//), [1](P,0,-I),C1]  A,(ti}-I)  ,C 

[7]  A+B+C+U+ ' ' 

[8]  DZ*-(lUHJETA)+((pZKlHnDBETA)o.+(pZKl'\*nDnETA+lffll 

V 
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uowcontingency 

ANALYSIS  OF  CONTINGENCY  TABLES 

(1)  CONTINGENCY  bi 

(2)  FOURFOLD  M 

(3)  MULTIPOLY  M 

(4)  Y+V  POOL  X 

THE  FIRST  TWO  FUNCTIONS  APPLY  TO  2-DIMENSIONAL  CONTINGENCY  TABLES ,  THE  OTHERS 
TO  CONTINGENCY  TABLES  IN  ANY  NU4BER  OF  DIMENSIONS.  ' CONTINGENCY '  PERFORilS  A  Cal 
SQUARED  TlST  OF  ASSOCIATION ,  WITH  DISPLAY  OF  STANDARDIZED  RESIDUALS.  ' FOURFOLD ' 
i4AY  BE  APPLIED  WHEN  THE  CATEGORIES  OF  EACH  CLASSIFICATION  ARE  ORDERED.  EMPIRI¬ 
CAL  LOG  CROSSPRODUCT  ( FOURFOLD )  RATIOS  ARE  DISPLAYED ;  A  PLACKETT  DISTRIBUTION  IS 
FITTED ,  AND  GOODNESS  OF  FIT  IS  TESTED  BY  CHI  SQUARED ,  WITH  DISPLAY  OF  STANDARDI¬ 
ZED  RESIDUALS.  'MULTIPOLY'  FINDS  EMPIRICAL  LOG  CROSSPRODUCT  RATIOS,  ANALOGOUS 
TO  THOSE  OF  ' FOURFOLD '  FOR  2 -DIMENSIONAL  TABLES.  'POOL'  POOLS  CATEGORIES  IN  ANY 
TABLE. 

(1)  CONTINGENCY .  THE  ARGUMENT  (M)  MUST  BE  A  MATRIX  OF  NONNEGATIVE  INTEGERS , 
HAVING  NO  ZERO  HARGInAL  TOTAL.  ( IF  M  HAS  FRACTIONAL  COUNTS ,  LINE  [1]  SHOULD 
BE  BYPASSED. )  GLOBAL  OUTPUT  VARIABLES:  EF  (MATRIX  OF  EXPECTED  FREQUENCIES),  SR 
(MATRIX  OF  STANDARDIZED  RESIDUALS).  SR  IS  DEFINED  THUS: 

SR  «-  (i-l-EF)*£F* 0.5 
+/,SR*SR  IS  CHI  SQUARED. 

NOTE.  FOR  GUTTMAN  PREDICTION  OF  THE  COLUMN  CATEGORY  FROM  THE  ROW  CATEGORY, 
LAMBDA  «-  ((+/r/M)-G)*(+/,W)-<M/+/A; 

(SEE  GOODMAN  AND  KRUSKAL:  JASA  49  (1954),  732-764.) 

(2)  FOURFOLD.  SAME  REQUIREMENT  FOR  M  AS  IN  'CONTINGENCY' .  THE  MATRIX  OF 
ICG  FOURFOLD  RATIOS  HAS  DIMENSION  VECTOR  (pfii)-l  ,  AND  REFERS  TO  ALL  POSSIBLE  2- 
BY-2  TABLES  THAT  CAN  BE  FORi-lED  FROM  M  BY  POOLING  ADJACENT  ROUS  AND  COLUMNS.  IF 

JiU.A.UU 

N 21  |  H22 

IS  SUCH  A  TABLE ,  THE  NATURAL- LOG  FOURFOLD  RATIO  (LFR)  CALCULATED  IS 
•+/0 . 5+0/11  ,N  12  ,N22  ,./21 ) 

THE  0.5  ADDED  TO  EACH  COUNT  HANES  LFR  A  ROUGHLY  UNBIASED  ESTItiATE  OF  LOG  PSI, 
WHERE  PSI  IS  THE  CORRESPONDING  FOURFOLD  RATIO  OF  PROBABILITIES,  AND  PREVENTS 
DISASTER  IF  AiJY  COUNT  IS  ZERO. 

POSSIBLY  THESE  VALUES  OF  LFR  DIFFER  ONLY  BY  SAMPLING  ERROR.  A  PLACKETT  DIS¬ 
TRIBUTION  tlAVING  CONSTANT  PSI  IS  FITTED  ITERATIVELY  BY  MAXIVM  CONDITIONAL  LIKE¬ 
LIHOOD ,  GIVEN  THAT  THE  MARGINAL  PROBABILITIES  ARE  PROPORTIONAL  TO  THE  GIVEN  MAR¬ 
GINAL  TOTALS.  AT  STAGE  0  THE  THREE  TRIAL  VALUES  FOR  LOG  PSI  ARE  THE  MOST  PRE¬ 
CISE  OF  THE  MEMBERS  OF  LFR  AND  THE  SAl-fE  ±  ONE  ESTIMATED  STANDARD  ERROR.  THE 
LOG  LIKELIHOOD  FUNCTION  IS  COMPUTED  FOR  THE  3  VALUES  OF  LOG  PSI  ,  A  PARABOLA  IS 
FITTED,  AiiD  AT  STAGE  1  THE  NEXT  3  TRIAL  VALUES  FOR  LOG  PSI  ARE  THE  VALUE  TO  MAX¬ 
IMIZE  THE  PARABOLA  AND  A  VALUE  ON  EITHER  SIDE,  USUALLY  CLOSER  THAI'  BEFORE.  AF¬ 
TER  AT  pIOST  4  ITERATIONS ,  EXPECTED  FREQUENCIES  (EF)  ARE  CALCULATED,  AND  HENCE 
STANDARDIZED  RESIDUALS  (SR)  AND  CHI  SQUARED. 

THE  4  ITERATIONS  WILL  NOT  FIND  THE  CORRECT  M.L.  VALUE  FOR  LOG  PSI  IF  THAT  IS 
INFINITE.  TO  ADJUST  THE  MAXIMUM  NUMBER  OF  ITERATIONS,  CHANGE  THE  '4'  IN  LINE 
[1]  TO  ANY  INTEGER  NOT  LESS  THAN  2. 

SEE 

PLACKETT:  JASA  60  (1965),  516-522. 

TUKEY:  EXPLORATORY  2dZd  /DIALYSIS.  ADD! SON  l/ESLLY ,  PRELIM.  ED.  1971,  CHAP.  29X. 
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(3)  MULiItuL,! .  TNE  hhiiUnLt.'i'  \i-j)  nJBT  nL  AH  ARRAY  UF  hum  LG  Ail ZL'  INTEGERS  IN 
2  OR  MORE  DIMENSIONS .  GLOBAL  OUTPUT  VARIABLES :  LCR  ( LOG  CROSSPRmCT  RATIOS), 
ESE  ( ESTIMATED  STANDARD  ERRORS).  LCR  AND  ESC  :!AVF  DIMENSTOH  VECTOR  (pM)-l  . 

IF  M  IS  A  MATRIX,  LCR  IS  THE  SAME  AS  LFR  YIELDED  BY  'FOURFOLD'.  IF  M 
IS  3-DIMENSIONAL,  LCR  REFERS  TO  ALL  POSSIBLE  2-BY-2-BY-2  TABLES  THAT  CM  BE 
FORMED  FROM  M  BY  POOLING  ADJACENT  PLANES,  HOI'S  AND  COLUMNS .  IF 

mi  i  uin  jmi-Liim _ 

N 121  |  i/122  A/221  |  i/222 

15  Si/67/  4  T4m',  THE  NATURAL-LOG  CROSSPRODUCT  ( EIGHTFOLD )  RATIO  LCR  IS 
»4/0.5+(i/lll  ,#112 , A/122, SI 21  ,N22\  ,.7211  ,/7212  ,N222 ) 

THE  VARIANCE  MAY  BE  ROUGHLY  ESTIMATED  BY 
+/f0.5+((/lll, i/112,  .  .  .  , i/222) 

TS£  SQ/MSS  SOOTS  ASS  Pi/T  Oi/T  /IS  SSS  ,  SIMILARLY  IF  M  HAS  MORE  THAN  THREE  DI¬ 
MENSIONS. 


(4)  POOL.  TNE  SECOND  ARGUMENT  ( X )  IS  AN  ARRAY  IN  2  OR  MORE  DIMENSIONS,  TYPIC¬ 
ALLY  A  CONTINGENCY  TABLE  OR  A  TABLE  OF  EXPECTED  FREQUENCIES.  TNE  FIRST  ARGUMENT 
(K)  IS  A  VECTOR  l/ITH  AT  LEAST  3  ELEMENTS.  r/[l3  SPECIFIES  THE  COORDINATE,  AND 
141/  SPECIFIES  THE  INDEX-VALUES ,  OVER  WHICH  THERE  IS  TO  BE  POOLING.  SECTIONS  OF 
X  CORRESPONDING  TO  INDEX-VALUES  2\V  OF  COORDINATE  Vlll  ARE  ADDED  TO  THE  SEC¬ 
TION  WITH  INDEX -VALUE  V[2  ]  ,  AND  THEN  THE  FORMER  SECTIONS  ARE  DELETED.  THUS  IF 
X  IS  A  MATRIX  WITH  5  ROWS,  AND  IF  V  IS  14  5  1,  POOLING  WILL  BE  OVER  THE  1ST 
COORDINATE  (ROWS);  THE  CONTENTS  OF  RO!'S  4,  5  AND  1  WILL  BE  ADDED,  PLACED  IN  ROW 
4,  AND  THEN  ROWS  5  AND  1  WILL  BE  DROPPED,  SO  THAT  THE  RESULT  (Y)  HAS  3  ROWS,  THE 
OLD  2ND  AND  3RD  AND  THEN  THE  SUM  OF  THE  OLD  1ST,  4 Til  AND  5TH  ROWS.  IF  V  HAD 
BEEN  114  5  THE  SAME  RESULT  WOULD  HAVE  BEEN  OBTAINED  EXCEPT  FOR  A  PERMUTATION 
OF  ROWS. 

9  JULY  1973 


V  CONTINGENCY  M;C;R 

[1]  *^2+  (  a/  ,  0  </•/ )  a  (  a  /  , M-  [/*/  )  a  (  a  /  2 <pA/ ) a  2 = p  p/  / 

[23  NO  GO,  THE  ARGUMENT  SHOULD  BE  A  MATRIX  OF  NONNEGATIVE  INTEGERS 

C  3  ]  *»4+  (  a  / 1  <0+  -f-M )  a  a  / 1  <R++ /M 

[4]  -p!>'SO  GO,  A  MARGINAL  TOTAL  VANISHES .  ' 

[5]  'THE  MATRIX  OF  EXPECTED  FREQUENCIES  IS  NAMED  EF 

[ 6 ]  ’ LEAST  MEMBER  OF  EF  -  ' ; L /  ,EF<-R °.*Cr+/S 

[  7  ]  '  STANDARD  TIED  RESIDUALS  (SR):' 

[83  0.01xL0.5+100xSS^(6'F*"0.5)xA/-SF 

[93  'CHI -SQUARED  =  '  ;+/,SSxSS;i>' ' 

[103  'DEGREES  OF  FREEDOM  =  ';x/  :+prf 
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7  FOURFOLD  M\U\R\C\GT\RCT\CCT\S\<J\D\L\ITS 
tlJ  -*,2+(  2£TTS«-4  )a(  a  /  ,A7>J«-0)a(  a/,/.;=|  A/)a(  A/2£pM)A2=ppM 

[2]  CO,  27/£  ARGUMENT  SHOULD  EE  /  HOI J NEGATIVE  INTEGER  MATRIX  WITH  POSIT 
IVE  MARGINAL  TOTALS. ' 

[3]  7M(itf)«>.>i/MpAf)[l])  +  .xA/+.x(i(>.Si(>(pAf)[2] 

[4]  GT+W[R;C] 

[5]  ^2xiv/(/fcr,cf)=o,/?cr^,((A?-i),'i)tw 

[6]  ->2xiV/(CC2,,C?,)=0,OJ2V,(~i,C-l)t!’ 

[7]  l/<-  1  "]  W 

[8]  ^.5+[/,((/?CTo.t(0-l)pO)-J/),(Gr+y-/?C2,°.+OOr),[2.5](((/?-l)pO)o.+CCr)-f/ 

[9]  'LOG  FOURFOLD  RATIOS  (LFR) : ' 

[10]  (L0.5+100x£Ftf-K.’T/,7)UOO 

[11]  LP+(tLFR)U/iS<-l/V+-,+/iW] 

[12]  LP+LP+  10  1  *S<-S+ 0.5 

[13]  kl-.'STAGE  LOG  PS  I  =  '  ;LP^-OP+lE"llxC!rxv/C!r>200000000000x  |LP 

[  14]  7+CT+  ( ~1+PSI+*LP)  o  .  x/fCTo  .  +OT 

[15]  fc4-(l/-(  (7x7)-(4xP5JxPSJ-l  )o  ,  xRCTa .  *CCT)*0. 5  )*(  2xPSJ-l )°  .  +  (pLFP)pO 

[16]  7-H7-0,[2]  0  “1  0  4-IM/-0,  0  0  "l  IIMf/,[2]( (pLP)pO)° .+CCT) , ((pLP)p0)°.+A?C2\ 
GT 

[17]  +[,3x\J=ITS 

[18]  £<-+/+/(  (pi/)pM)xffi,V 

[19]  ->£2xi02PW.L[2  12  3] 

[20]  Lr4.,LP[2>ZMS*-/L[3  l])*2xP 

[21]  +lilx\ITS=J+J+l 

[22]  I2SM1  0  =S22x|P)/(J+l),ir5 

[23]  -£,l,p  LWP+  "10  1  xS-H(CT4)f(0.5x|P)LS 

[24]  12:' LIKELIHOOD  FN  CURVES  WRONG  WAY.  ENTER  3  INCREASING  EQUAL-SPACED  VALUES 

FOR  LOG  PSI. ' 

[25]  +(&lxl((-/“2+LP)=-/2tLP)A0>-/2+LF4-3+,0),&2 

[26]  Hi:' THE  MATRIX  OF  EXPECTED  FREQUENCIES  IS  NAMED  EF  .  \0" 

[27]  'LEAST  MEMBER  OF  EF  =  ’  ;L/,EF«-f/.V 

[28]  ' STANDARDIZED  RESIDUALS  (SR):' 

[29]  0.01x|_0.  5+100xS/?+-(E'F*-0. 5)*M-EF 

[30]  'CHI -SQUARED  -  ' ;+/ ,SR*SR ; i>" 

[31]  'DEGREES  OF  FREEDOM  =  ';"l+x/p LFR 

V 


7  MULTIPOLY  MiliKy’ 

[1]  *(pl-*-  1  0 )+  ( A/  ,/./s0  )  A  (  a/  ,/;  =  LM )  A  (  a  /  2  Sp/'  / )  A  2  <K<-p  pM 

[2]  -K),p Q*-'NO  GO.' 

[3]  -*-3x\(2*/0>pI^I,~P 

[4]  -+4x  i(2x/0>PPW-«-(^ ;  (opM)-K) ,  ((ppM)-K)+(2lr\K+l ) ,  2  l)^+.x(~P),[ 
2.5]  //«-( i//)» .  >i  l+A'-*-"itp// 

[  5  ]  AM ,V«-1 ) + 2>  ( (  pi ) , /a pM )  p/f 

[6]  tfj'-WLPimF  INTERACTION— LOG  CROSSPRODUCT  RATIOS  (LCR) : ' 

[7]  0 . 01  x  L  0 . 5+1 00  xLCP-?-o(  x  /  J /A/ )  *  x  /  (~  j )  /<.; 

[8]  'ESTIMATED  STANDARD  ERRORS  (ESE):'  ,[>' ' 

[9]  0.01xL0.5+100xPS^-(?x  +  /fA/)*0.5 
V 
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CHARACTERISTIC  ROOTS  AND  VECTORS  OF  A  SYMMETRIC  MATRIX 
Y*C  JACOBI  X 

THE  SECOND  ARGUMENT  (X)  IS  THE  GIVEN  MATRIX,  TO  BE  TRANSFORMED  TOl/ARDS  DIAGO¬ 
NAL  FORM  BY  JACOBI'S  METHOD.  THE  FIRST  ARGUMENT  (C),  A  POSITIVE  SCALAR.  IS  THE 
TOLERANCE  FOR  OF I  DIAGONAL  ELEMENTS  IN  THE  TRANSFORMED  MATRIX.  THE  EXPLICIT  RE¬ 
SULT  (Y)  IS  A  THREE-DIMENSIONAL  ARRAY  IfITH  2  PLANES.  Y[l;;]  IS  THE  TRANSFORMED 
MATRIX,  AND  SO  1  1  6)Y[l;;]  IS  THE  VECTOR  OF  ROOTS.  UHICll  IS  ALSO  OUTPUT  AS  THE 
GLOBAL  VARIABLE  'ROOTS'.  THE  CORRESPONDING  VECTORS  ARE  THE  COLUMNS  OF  Y[2;;]  . 
THE  GLOBAL  SCALAR  VARIABLE  'ITERATIONS'  IS  THE  NUMBER  OF  PERFORMED  STEPS  OF  THE 
JACOBI  PROCESS. 


A 


USUALLY  (IF  THERE  HAVE  BEEN  ENOUGH  ITERATIONS)  THE  ROOTS  APPEAR  IN  DECREASING 
ORDER. 

TO  SEE  THE  EFFECT  OF  REDUCING  THE  VALUE  OF  C  ,  LET  Cl  AND  C2  BE  TOLERAN¬ 
CES,  THE  FIRST  GREATER  THAN  THE  SECOND.  ENTER.  FOR  A  GIVEN  MATRIX  X  , 

Y1*C1  JACOBI  X 
ROOTS 

ITERATIONS 

Y2+C2  JACOBI  YlCl ; ; 3 
ROOTS 

ITERATIONS 

THE  SECOND  VERSION  OF  'ROOTS'  IS  THE  BETTER  ONE.  THE  CORRESPONDING  IMPROVED 
EIGENVECTORS  ARE  THE  COLUMNS  OF  Yl[2 ;;  3+ .  *Y2[2 ; ;  ]  . 

5  DEC.  1972 


V  Y+C  JACOBI  X\A-,I\J-,N\M-,N-,R-,S-,Z 

[1]  ^2+(A/C>0)A(0=ppC’)A((t/,C)>r/|(,Y)-,«pX)A(4<+/pX)A(=/pY)A2=ppY 

[2]  '♦O.pO-'/VO  GO.  '  ,Y«-'  ' 

[3]  !«-,(  \N)°  .  <\N+-+/ltpY+-X 

[4]  R+(N,N) pi, NpITERATIONS+O 
[53  Li:-'(Ca‘M/Z+|//,Y)/£4 

[63  J+U-(M-K+-l+ll\  l+M«-(  J\Z ) \M)*N 

[73  -(Y[e/;/3*Y[K;A3>/£.2 

[83  ■+L3,p/l-*-o  0.25xxY[J;/G 

[93  £-2:/l<-0. 5"(or[«7;J3<y[A;tf3)+  302*YW;A'3*Y[A;A'3-Y[<7;  </3 

[103  43:S«-(tf,«)pl,iVp0 

[113  SLK-,K,jy-  1  *1  *5lJ-,J,K)<-  2  1  °.oA 

[123  YH$R )+ .  x.Y- .  */,<-/?+ .  xs 

[1 3  3  +ftl ,  p/r;  71 1  IONS*- ITER  A I  IONS + 1 

[143  £4  -.ROOTS*-  1  1  i?Y 

[153  Y<-Y,[0.53  R 
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HOWPLOT 

THREE  FUUCTIOilS  FOR  PLOTTING 

(1)  IMI  SCATTERPLOT  V 

(2)  DOUNPLOT  V 

(3)  S  RGWPLOT  V 

'SCATTERPWT'  IS  USED  FOR  JCATTERPLOTTING  CORRESPONDING  MEt IDLES  OF  TWO  VEC¬ 
TORS.  THE  WHOLE  CHARACTER  ARRAY  IS  GENERATED  IN  THE  WS  AND  MAY  THEN  BE  PRINTED. 
'DOUNPLOT'  IS  USED  FOR  PLOTTING  MEMBERS  OF  ONE  OR  MORE  VECTORS  AGAINST  THEIR  IN¬ 
DEX  NUMBERS.  THE  CHARACTER  ARRAY  IS  GENERATED  AiJD  PRINTED  LINE  BY  LINE ;  IT  IS 
NOT  STORED  IN  THE  WS.  'ROW PLOT'  GENERATES  AtID  DISPLAYS  IN  ONE  ROW  THE  CONTENTS 
OF  EACH  OF  ONE  OR  MORE  VECTORS.  SEE  ALSO  ' UOUTRIPLEPLOT' . 

(1)  SCATTER -LQT.  U  AND  V  ARE  VECTORS  OF  EQUAL  LENGTH.  CORRESPONDING  MEMBERS 
UIJZ  AND  Vl-JZ  /HE  PLOTTED  AS  ABSCISSA  AND  ORDINATE  OF  A  POINT.  A  SINGLE  POINT 
IS  SHOWN  AS  o,  TWO  COINCIDENT  POINTS  AS  9,  THREE  OR  MORE  COIHCIDENT  POINTS  AS  •. 

( YOU  MAY  CHANCE  THIS  CODE  E.  JSILY — SEE  BELOW. ) 

TO  GE *  A  SINGLE  PLOT  DO  NOT  NAME  AN  EXPLICIT  RESULT  BUT  CALL  THUS : 

U  SCATTERPLOT  V 

TO  GET  A  PLOT  AND  ALSO  STORE  THE  SYMBOL  ARRAY  FOR  MAKING  A  COPY  LATER  CALL  THUS: 
Of hV  SCATTERPLOT  V 

THE  USER  IS  INTERROGATED  ABOUT  THE  SIZE  OF  THE  PLOT  (NO.  OF  ABSCISSA  VALUES, 
NO.  OF  ORDTir.TE  VALUES ) ,  AND  ABOUT  WHETHER  HE  WANTS  THE  HORIZONTAL  AND  VERTICAL 
SCALES  Tu  BE  EQUAL.  IF  THE  SCALES  ARE  NOT  MADE  EQUAL  THEY  ARE  CHOSEi  LIKE  THIS. 
THE  ABSCISSA  Vj  RUE  SET  IS  EQUAL- SPACED ,  THE  LEAST  EQUAL  TO  (l /{/),  THE  GREATEST 
TO  (1710.  EXCEPT  THAT  IF  THESE  EXTREMES  DO  NOT  HAVE  THE  SAME  SIGN  THE  VALUE  SET 
IS  MINIMALLY  TRAilSLATED  WITHOUT  CHANGE  OF  SCALE  SO  THAT  0  BECOMES  A  MEMBER  OF 
THE  SET.  SIMILARLY  FOR  ORDINATES.  IF  THE  SCALES  ARE  MADE  EQUAL.  EITHER  THE  AB¬ 
SCISSA  VALUE  SET  OR  THE  ORDINATE  VALUE  SET  MAY  BE  SHORTENED,  TO  AVOID  UNNECESSA¬ 
RY  BLANK  COT  WINS  OR  ROWS  IN  THE  PLOT.  THE  COORDINATES  OF  THE  PLOTTED  POINTS  ARE 
ROUNDED  U  Tir  NEAREST  MEMBERS  OF  THE  COORDINATE  VALUE  SETS.  THE  AXES  ARE  PRIN¬ 
TED  AS  *'S.  EXCEPT  THAT  IF  0  OCCURS  IT  IS  SO  SHOWN. 

TO  CHAtrC  Tn „  PLOTTING  CODE  ALTER  THE  CHARACTER  VECTOR  IN  LINE  Cl].  FOR  EX¬ 
AMPLE,  TO  L'iOU  ITJLTIPLI CITIES  BY  NUMERALS,  1  TO  9.  WITH  X  FOR  10  OR  MORE,  RE¬ 
PLACE  'eeo'  NY  '*0870 5“ 321'.  (MC  SHOULL  NOT  CONTAIN  A  BLANK  OR  \  OR  -.) 

THE  TERMINAL  TS  SUPPOSED  SET  *T  SINGLE  SPACING .  AlID  PAGE  WIDTH  AT  THE  MAXIMUM 
OF  130  CtM  CTENJ  IF  NEEDED.  THE  LAST  LINE  (2  A  SPECIFIES  A  HORIZONTAL  SPACING 
OF  THE  DIi  7/7;  YOU  M\Y  DELETE  THIS  LINE ,  IN  WHICH  CASE  THE  *64'  IN  LINE  [7] 
SHOULD  BE  C'AVED  TO  '129'.  THE  LOCAL  VARIABLE  S  IN  LINE  [12]  SHOULD  GIVE  THE 
RELATIVE  MIEN  ETUDES  OF  THE  HORIZONTAL  AND  THE  VERTICAL  UNIT  STEPS,  AS  PRINTED  BY 
YOUR  TER  INAL.  WITH  6 -PER- INCH  LINEFEED,  THE  CORRECT  SETTINGS  FOR  S  ARE 
:on  12-PITCH  CHARACTER  SPACING,  LAST  LINE  INTACT:  S+-  1  1 

l OR  10-PITCH  CHARACTER  SPACING,  LIST  LINE  INTACT:  S+-  6  5 

IIP  ±2 -PITCH  CHARACTER  SPACING,  L/LT  LINE  DELETED:  S«-  12 

FOR  10 -PITCH  CHARACTER  SPACING,  LAST  LI  IS  DELETED:  S«-  3  5 

(SAME  FOR  'T3CP'  ,  THE  LAST  LINE  BEING  NUMBERED  [33].) 

(2)  do': ’/  is  either  (i)  a  vector  or  (2)  a  matrix,  if  (1),  viJZ  is 

PLOTTED  AGAINST  .1 .  IF  (2).  THE  ROWS  ARE  PLOTTED,  I.E.  VlM  AGAINST  J.  POINTS 
ARE  SHOWN  Ro  MULTIPLICITY  iJ  NOT  INDICATED.  •*-  OR  -*  MEANS  AT  LEAST  ONE  VALUE 
WAS  OUTER DM  THE  PLOTTING  RANGE. 


THE  USER  .  "jHT  SPECIFY  SCALE  AND  LABELING  (WITH  THE  GLOBAL  VARIABLES  DPS  AND 
DPL  ),  AMD  IS  OFFERED  INSTRUCTIONS. 
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TO  CHMGE  THE  SYMBOLS  (o,  -)  USED  IN  PLOTTING ,  ALTER  LINES  [33]  AND  [37]. 

21?  CHANGE  THE  TOP  AND  BOTTOM  ROUS  OF  THE  PLOT,  14ARKING  THE  ABSCISSA  SCALE ,  ALTFR 
THE  DEFINITION  OF  AC  IN  LINE  [27].  LABELING  OF  THE  TOP  ROM  CAN  BE  INSERTED  AS 
A  NEW  LINE  [26.1]. 

(3)  ROUP  LOT.  S  IS  A  VECTOR  SPECIFYING  SCALE ,  JUST  LIKE  DPSi  SEE  INSTRUCTIONS 
IN  'DOWNPLOT' .  V  IS  EITHER  (1)  A  VECTOR  OR  (2)  A  MATRIX.  IF  (1).  THE  SET  OF 
VALUES  IS  SHOWN  AS  POINTS  ON  A  LINE.  IF  (2),  EACH  ROM  IS  SO  SHOWN ,  ONE  LINE  PER 
ROW.  MULTIPLICITY  IS  SHOWN  AS  IN  ' SCATTERPLOT ' .  TO  CHANGE  THE  PLOTTING  CODE 
ALTER  THE  CHARACTER  VECTOR  IN  LINE  [4]. 

25  FEB.  1972 


V  W+U  SCATTERPLOT  V\E\I',J -,M\MC\N',R-,S 

[1]  -*-2+(i=ppM>,«eo'  )A(i=ppU)A<i=ppi')A(i(te2)A(*/py)=//*-+/py 

[2]  -+pO- ' ARGUMENTS  SHOULD  BE  VECTORS  OF  EQUAL  LENGTH  NOT  LESS  THAN  2 

[3]  -4+A/0 <R*~/E+  2  2  /U)Al/U),(\ /V),UV 

[4]  -HJ.pO  W  ARGUMENT  HAS  ZERO  RANGE.  '  ,V<-' ' 

[5]  'SIZE?  TYPE  TWO  NUMBERS,  SUCH  AS:  25  25' 

[6]  +5*\(v/M*\.M)v(v/lZM)v2*pM+,L\ 

[7]  **8+64£A/[l] 

[8]  -*•6,p[>,FJ/^S2,  NUMBER  TOO  LARGE,  TRY  AGAIN.' 

[9]  I+R+M-J+l 

[10]  'SAME  SCALES?  Y  OR  N.  ' 

[11]  -K(,/VY'=lt,H])/  13  12). 10 

[12]  £[;2>E[;l]-I*“UAM+r/?*I[l  2  ]<*£*[//♦$«•  1  1 

[13]  (0W) p*  ').[1] 

[14]  +£1M0<*/  1  2  t£4*F*2’,[l. 5]  I 

[15]  E[l;2Xl-M[l])+E[l;lM0.5+£[l;l] 

[16]  I/[;(pf/)[2j-ff[l;l]MM[2]p'|').’0' 

[17]  £1  'EXTREME  ABSCISSAS  ARE:  '  ;$I[1>£[1;] 

[18]  'ABSCISSA  UNIT  STEP  IS  '  ;ltl 

[19]  ->£2M0<x/£[2;] 

[20]  F[2;2Wl-W[2])+£,[2;l>L0.5+ii’[2;l] 

[21]  */[l+ff[2;l];>»0l,*[l]p'-' 

[22]  U2 :  'EXTREME  ORDINATES  ARE:  '  ;<H[2]*ZT[2;] 

[ 23 ]  ' ORDINATE  UNIT  STEP  IS  '  ;I[2];£,£ 

[24]  Xp4')[2]-rF[l;l]-0.5+i/*I[l] 

[25]  V*\  0 . 5+fit  2 ;  1]  -  *I[  2  ] 

[26]  £3:(/[l'[J];i/[J]XC[ir~l+WCi*V[l'[«/];i/[^]]] 

[27]  +L3*\N*J<-J+l 

[28]  IPrV*" 

[29]  M+(  ( 2x_ltpJ/)p  0  1)W 

V 
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7  DOWNPLOT  V-,J-,N-,C-,D-,AC 
Cl]  -*■((( A/lSplf)A  1  2  =ppV)/  3  4), 2 

C2]  -*0  ,pQ*-'  ARGUMENT  WON"?  DO.  ' 

C3]  Nl.pW 

[4]  'INSTRUCTIONS?  V  OR  IV 

[5]  -*-((’i/r=lt,U)/  21  6), 4 

[6]  'TWO  GLOBAL  VARIABLES  MUST  BE  SPECIFIED,  DPS  ( DOWNPLOT  SCALE)  ADD  DPL  ( D 
OWNPLOT ' 

[7]  'LABELING).  BOTB  ARE  VECTORS.  DPS  INDICATES  THE  PLOTTING  SCALE  FOR  ARGUM 
BUT  VA-X 

[8]  'LUES,  LISTING  THE  LOWEST  VALUE ,  THE  STEP  SIZE  AND  THE  HIGHEST  VALUE  TO  BE 
ACCOM-' 

[9]  'MODATED.  HIGHEST  MINUS  LOWEST  MUST  BE  DIVISIBLE  BY  STEP.  EVERY  TENTH  V 
ALUS  IS' 

CIO]  'MARKED  BY  +,  UNLESS  LINE  C27]  IS  ALTERED.  ENTER  E.G.  (45  PLOTTING  VALUES 
):' 

[11]  •  DPS*-  8.1  0.1  12.5’ 

Cl 2]  'DPL  INDICATES  THE  LABELING  OF  THE  INDEX  NUMBERS,  LISTING  FOUR  INTEGERS. 
THE  1ST' 

C 1 3 ]  'MUST  BE  POSITIVE  AND  NOT  GREATER  THAN  THE  3RD.  THE  DPLllV'TH  ROW  OF  THE 
PLOT  IS' 

C14]  'LABELED  DPL[2) ,  AND  THEREAFTER  EVERY  DPLl3V'TH  ROW  IS  LABELED  WITH  AN  I 
NCREMENT' 

Cl 5]  'OF  DPLM.  FOR  SIMPLE  COUNTING  OF  ROWS  BY  TENS  ENTER-.' 

Cl 6]  '  PPL- 4pl0‘ 

C17]  'WITH  DPS  AND  DPL  SPECIFIED,  RESUME  EXECUTION  OF  DOWNPLOT  BY  TYPING-.' 

Cl8]  '  ->20' 

[19]  SADOWNPLOT*~ 9+ p|> » SEff  HOWPLOT  FOR  FURTHER  INFO.  ' 

[20]  SADOUNPLOT- K) 

[21]  +23-^v/(3=pDPS)A(A/O*~im,0PS)Al=pp/?PS 

[22]  -►23+(//=U)a(125^)a2S/^1+(-/PPSC3  l])+PPS[2] 

[23]  -0,p|> 'DPS  WON'  "T  DO.  ' 

[24]  -*26-(A/,£PL=l0PP)A(4=+/pPPZ,)Al=pp£>PP 

[25]  ->26+(UPPL[1])a>/PPL[3  1],0" 

[26]  -O.pO'PPL  WOtV'T  DO.' 

[27]  »  +'  AN+hDpAC*-'- . +  ' 

[28 j  £1  :C-(/7+2)pl 

[29]  C[ir(L2.5+(7[;«/]-DPS[l])*PPS[2])L//+2]-2 

[30]  C[2]-3LC[2]+2*CT  1  ]-l 

[31]  C[tf+l]-4lC[/i/+l]+3xC’U'+2]-l 

[32]  -&3xi0=PPL[3]|J-PPL[l] 

[33]  '  I','  O—'tU'UC] 

[34]  £2:-&l*i(py)[2]>J-J+l 

[35]  —0,pL}-'  +'  ,NpAC 

[36]  L3:ZM4pl0)TPPL[2]+PPL[4]x(J-DP£[l])*PPL[3] 

[37]  -£2,pO-'0123456789  ' [l+P+(  10x(P=0)a.  v( \4 )o .  >i3) ,0] , '  |»  , ' 
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HOW I  INTEGRATE 

ONE- DIMENSIONAL  DEFINITE  INTEGRALS 
ZHI  INTEGRATE  A 

THE  FIRST  ARGUMENT  (//)  IS  A  STEP  SIZE,  SCALAR.  THE  SECOND  ARGUMENT  (/l)  IS  A 
VECTOR  OF  2  OR  MORE  LIMITS  OF  INTEGRATION,  IN  ASCENDING  ORDER,  WITH  DIFFERENCES 
ALL  DIVISIBLE  BY  H  .  THE  EXPLICIT  RESULT  (Z)  IS  A  VECTOR  OF  LENGTH  1  LESS  THAN 
THE  LENGTH  OF  A  ,  LISTING  THE  DEFINITE  INTEGRALS  FROM  /i[l]  TO  EACH  OF  THE  OTHER 
MEMBERS  OF  h  .  THE  FUNCTION  TO  BE  INTEGRATED  IS  ASKED  FOR,  AND  MUST  BE  EXPRES¬ 
SED  IN  TERMS  OF  AN  ARGUMENT  X  ( LOCAL  VARIABLE).  TWO-POINT  GAUSSIAN  QUADRATURE 
IS  USED  IN  EACH  INTERVAL  OF  WIDTH  li  . 

FOR  EXAMPLE ,  TO  INTEGRATE  SIN  X  FROM  0  TO  EACH  OF  12  3  AND  COMPARE  THE 
RESULT  WITH  1  -  COS  X  : 

0.1  INTEGRATE  “l+i4 
FUNCTION  JF  X  TO  BE  INTEGRATED? 

J: 

lox 

0.4596976835  1.416146804  1.909992451 

1-20x3 

0.4596976941  1.416146G37  1.989992497 

29  NOV.  1972 


V  Z +11  INTEGRATE  A  \N\X 

[1]  -*-3-(A/2<Pi4)A(l=pp/l)A(A/,/  >0)A0=Pp// 

[2]  ->3+(A/0<i7-0,“l+^)AA//;=L4'-<-((  UA  )-!+/*  )m 

[3]  ->0,p[>'/W  GO.  '  ,Z*-'  ' 

[4]  ' FUNCTION  OF  X  TO  BE  INTEGRATED ?' 

[5]  tf*-UCl ]+//* ri+/l/)o.-('/i2)>l+  1  "l  *3*0.5 

[6]  Z<-UliA)*(N°.>\'l  M)+.x+/L 

V 


V  S  ROWPLOT  V\J  \N  ■,(>,!  1C 

[i]  ^(((A/l<pI/)A  1  2  =ppl7)/  3  4  )  ,  2  0 

[  2  ]  -pt> '  SECOND  ARGUMENT  WON  "T  DO .  ' 

[3]  V-(l,pV)pV 

[  4  ]  -»6  • -  ■ v  /  ( 3 = p/; )  a  ( l = p  p5 )  A 1 = p  pM7«|>  '  *90 ' 

[ 5 J  +6+(125^  U/W//-LJ>A2<iKl  +  (-/5[3  l])*S[2] 

[  6  ]  +pi> '  FIRST  Ah  GUNS' ’7  WON  "T  DO.' 

[7]  C-*-(p.’:;C)L  +  /(  . =11  ( 12. 5+(  VI J ;  ]-5[l ] )*S[2]  )l//+2 

[8]  C[ 2 ]■*-(()  1  =1[HC)/C[2] ,  2 

[9]  CC/'+lMO  1  =  ll-ltC)/CU7+l],-l 

[10]  5p 1  V/CK3m  1+C] 

[11]  -*-7x  i  ( ltpF/'J-'  ,/+l 

[12]  5p *  !  ;//p  1 - +' 

V 
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FUNCTIONS  FOE  PLOTTING 

(1)  Z  TDP  V 

(2)  U*V  TSCP  V 

(3)  ROWCOLDI SPLAY  I 

IN  THESE  VERSIONS  OF  FUNCTIONS  DESCRIBED  Ili  'HOUPLOT' ,  A  THIRD  DIMENSION  IS 
SUGGESTED  BY  VARYING  THE  SYMBOL  USED  IN  PLOTTING  THE  POINTS. 

(1)  TQE  (TRIPLE  DOWNPLOT)  IS  A  MODIFICATION  OF  'DOUNPLOT'.  THE  GLOBAL  VARIA¬ 
BLES  DPS  AND  DPI  MUST  BE  SPECIFIED  BEFORE  EXECUTION.  THE  FIRST  ARGUMENT  (Z) 
IS  A  CHARACTER  SCALAR  OR  1 VECTOR  SHOOING  WHAT  SYMBOLS  ARE  TO  BE  PLOTTED.  IF  Z  IS 
SCALAR,  THE  SAME  SYMBOL  WILL  BE  PLOTTED  EVERY  TIME.  THE  CALL 

»o*  TDE  V 

YIELDS  WE  SAME  RESULT  AS  THE  CALL 
DUNNE  LOT  V 

WITH  THE  EXCEPTIONS  THAT  INSTRUCTIONS  ARE  NOT  OFFERED  AND  THAT  IF  ANY  MEMBERS  OF 
V  FALL  OUTSIDE  THE  RANGE  SPECIFIED  BY  DPS  THEY  ARE  MISSED  (NO  WARNING  OR  -*•). 

IF  Z  IS  NO?  SCALAR,  IT  MUST  BE  A  CHARACTER  VECTOR  OF  LENGTH  EQUAL  TO  ~UpV  . 
ZC«H  IS  THE  SYMBOL  USED  IN  PLOTTING  V£«7]  (IF  V  IS  A  VECTOR)  OR  Vt  ;«H  (IF  V  IS  A 
MATRIX)  AGAINST  J . 

TO  CHANGE  THE  TOP  AND  BOTTOM  ROWS  OF  THE  PLOT  i'iARKING  THE  ABSCISSA  SCALE,  AL¬ 
TER  THE  DEFINITION  OF  AC  IN  LINE  [12].  LABELING  OF  THE  ‘TOP  ROW  CAN  BE  INSERT¬ 
ED  AS  A  NEW  LINE  [11.1]. 

(2)  TSCP  ( i RIFLE -SCATTERPLOT )  IS  A  MODIFICATION  OF  'SCATTERPLOT' .  THE  ARGU¬ 
MENTS  (U,  V)  SFCCIFY  ABSCISSAS  Ai.’D  ORDINATES,  AS  BEFORE.  ALL  COINCIDENT  POINTS 
ARE  SHOWN  AS  *  .  AFTER  THE  INTERROGATIONS  ADOUT  FINENESS  OF  PLOTTING  AND  EQUAL 
SCALES,  TiIE  USER  IS  ASKED  FOR  THE  SYMBOLS  TO  BE  USED.  IF  HE  REPLIES  WITH  OHE 
SYMBOL  IN  QUOTES,  SUCH  AS  'o' ,  THE  RESULT  WILL  BE  LIKE  THAT  OF  ' SCATTERPLOT '  EX¬ 
CEPT  THAT  ALL  MULTIPLE  POINTS  ARE  SHOWN  BY  THE  ONE  SYtSOL  *  .  OTHERWISE  HE  MUST 
REPLY  WITH  A  CHARACTER  VECTOR  OF  LENGTH  EQUAL  TO  p U  AND  pV  ,  NONE  OF  THE  CHA¬ 
RACTERS  BEING  A  HI  AUK  OR  *  .  FOR  A  VECTOR  V  ,  THE  CALL 

(ip7)  TSCP  V 

WITH  Z  AH  ?  r  AUSWSU  TO  THE  THIRD  QUESTION  YIELDS  A  PLOT  EQUIVALENT  TO 
Z  TDP  V 

IF  THE  SCALES  ACE  SUITABLY  HATCHED.  BECAUSE  THE  MINUS  SIGN  MAY  BE  NEEDED  AS  A 
PLOTTING  SYMBOL,  ZEEO  HUES  ARE  NOT  INDICATED  IN  THE  PLOT  WITH  |  OR  -  ,  AS  THEY 
ARE  IN  '  SCATTERPLOT  • . 

(3)  &,iU7':>.TS?J.AK  IS  A  SPECIAL  FUNCTION,  SUBSTITUTING  FOR  'TSCP' ,  THAT  MAY  BE 

USED  TO  RILTIA.Y  THE  OUTPUT  OF  'ROW COL  ’ .  ABSCISSAS  AND  ORDINATES  ARE  THE  COLUMN 
EFFECTS  ATM  THE  R'-W  EiFECTS.  THE  SYMBOLS  PLOTTED  REPRESENT  THE  RESIDUALS  GRADED 
ON  A  1 -POINT  SCALE  FRCM  URGE-NEGATIVE  TO  LARGE-POSITIVE:  P,  £. 

AT  THE  H.yPFHSE  OH  ACCURACY  IN  POSITIONING  THE  PLOTTED  POINTS,  COINCIDENCES  OF 
ROWS  OR  COLUM  S  ATE  AVOIDED  THROUGH  A  UNIT  MINIMUM  SEPARATION,  EVEN  WEN  TWO  ROW 
EFFECTS  OR  COLUMN  EFFECT'S  ARE  EQUAL. 

THE  ARGUMENT  (I)  IS  THE  CHANGE  IN  COLUMN  EFFECT  REPRESENTED  BY  A  UNIT  HORI¬ 
ZONTAL  DISPLACEMENT.  TRY  EXECUTING  THE  FUNCTION  WITH  A  LARGISH  VALUE  FOR  I,  NOT 
LESS  TUAN  ((i /CF)-l/CE)l 20  .  IF  THE  DISPLAY  IS  TOO  SMALL  REPEAT  WITN  A  SMALLER 
VALUE  FOR  I 
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THE  FUNCTION  REFERS  TO  THE  GLOBAL  VARIABLES  CE,  RE,  RY  GENERATED  BY  'RO'/CUL', 
BUT  THESE  SHOULD  FIRST  BE  REARRANGED  BY  CALLING  THE  FUNCTION  'ROUCOLPERMUTE'  (HO 
ARGUMENTS  OR  EXPLICIT  RESULT),  WHICH  PERMUTES  THE  RONS  AND  COLUMNS  SO  THAT  CE  IS 
IN  ASCENDING  ORDER  AND  RE  IN  DESCENDING  ORDER. 

ADJUSTMENTS  III  LINE  Cl]:  TO  SUPPRESS  HORIZONTAL  DOUBLE  SPACING  SET  (D+l). 
FOR  b- PER -INCH  LINEFEED,  12 -PITCH  CHARACTER  SPACING,  IF  D+2  THEN  SM  1  ,  BUT 
IF  O+l  THEN  >12,  FOR  IQ -PITCH  SPACING,  SET  S+  6  5  AND  >35,  RESP. 

2 P  FEB.  1972 


V  W+U  TSCP  V\E\I\J\M\N\R\S\^ 

[1  ]  •*,2+(  l=pp£/)A(  l=ppl/)A(/7>2  )a  ( +/p£)  )=/i/«-+  /  qV 

[2]  +pD+' ARGUMENTS  SHOULD  BE  VECTORS  OF  EQUAL  LENGTH  NOT  LESS  THAN  2.  ' ' 

C3]  -4+a/O <R+-/E+  2  2  p(t/U),(l/U),(r/V),\./V 

[4]  -K),pl R-'ONE  ARGUMENT  HAS  ZERO  RANGE.  '  ,W+' ' 

[5]  'SIZE?  TYPE  TWO  NUMBERS,  SUCH  AS:  25  25' 

[  6  ]  -*•  5  x  i  ( v  /M  *  L  M )  v  ( v  / 1  HI )  v  2  *  &  -1+ ,  □ 

[7]  -8+64^[l] 

[8]  +b,p[}+' FIRST  NUMBER  TOO  LARGE,  TRY  AGAIN.' 

[9]  I+RrM-J+1 

[10]  'SAME  SCALES?  YOU  11.' 

Cll]  -*>((  Wslt.C])/  13  12),  10 

[12]  ^;2>tf[;l]-Ix~l+AM+r/?*I[l  2>5x[/Jt5-  1  1 

[13]  fc'-'x\((<tW)p'  «),[1]  ’x* 

[14]  ->^lxiO<x/  l  2  tE+E*I,L  1.5]  I 

[15]  ai;2Hl-wtl]Ml;lKo.5+C[l;l] 

[16]  </[  1 +p{/;  ( p//)  [  2  ]  -i,’[  1 ;  1 3  ]^ '  0 ' 

[17]  LI:' SYMBOLS  TO  BE  USED?' 

[18]  -*-42xi0<x/£[2;] 

[19]  £[  2 ;  2  ]•♦*(  1  -M[  2  ]  )+£[  2 ;  1 J+-L  0 . 5+£[  2 ;  1  ] 

[20]  (/[.lt/?[2  ;  1]  jlD-i-'O  ' 

[21]  42  :-*•(( (~v/,Zc'  *')'v(a///=p£)a  0  1  =pp4^D)/43+  1  2 ) ,43 

[22]  43:-*-(43-l),p[>'/W  GOOD,  TRY  AGAIN.' 

[23]  Z+NpZ. 

[ 24 ]  C,' EXTREME  ABSCISSAS  ARE :  '  ;4>I[  1  ] x£[  1 ;  ] 

[25]  'ABSCISSA  UNIT  STEP  IS  ' -,UI 

[26]  'EXTREME  ORDINATES  ARE:  '  ;<t >I[  2  ]x£[2 ;  ] 

[27]  'ORDINATE  UNIT  STEP  IS  '  ;I[2];G,£ 

[28]  U+( pW)[2]-ffc’[l : l]-0. 5+4*7[l] 

[29]  V+f  0.5tE[2;l]-I7*J[2] 

[30]  [A-MVtJhUlJWHWlvlJhUlJ)-)*'  ')/'*'),ZUl 

[31]  ■+44x  iN^J+J+l 

[32]  U+V+Z.+  ' ' 

[33 J  f/*-( ( 2x~ltp(/)p  0  1)W 

V 
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uoimx 

maximization  of  a  function  of  one  variable 
z+x  max  y 

THE  ARGUMENTS  ARE  VECTORS  OF  LENGTH  3,  THE  FIRST  HAVInG  NO  TWO  MEMBERS  EQUAL. 
THE  EXPLICIT  RESULT  IS  THE  COORDINATES  OF  THE  VERTEX  OF  A  PARABOLA  WITH  VERTICAL 
AXIS,  THAT  GOES  'THROUGH  THE  3  POINTS  WHOSE  ABSCISSAS  ARE  THE  FIRST  ARGUMENT  AND 
ORDINATES  THE  SECOND  ARGUMENT.  IF  THE  PARABOLA  IS  CONVEX  FROM  BELOW,  THE  VERTEX 
IS  THE  MINIMUM  INSTEAD  OF  THE  MAXIMUM  POINT,  AND  THE  WORD  'MINIMUM.  '  APPEARS  IN 
THE  OUTPUT. 

20  NOV.  1972 


V  Z+X  MAX  Y\A\B’,M 

Cl]  -3-(3rt/pX)A(3=+/py)A(i=ppX)Ai=ppy 
[2]  -*-3t(  V(H’/)=^U)a=/2^  0  0 

C  3  ]  -»pO  ’  ARGUMENTS  WON  "T  DO.'  ,Z+" 

[4]  >MflSl,X,[1.5]  X*X+X-tM+/X)i 3 

[5]  -»(U[3]=0),/i[3]<0)/  9  7 

[6]  'MINIMUM. ' 

[7]  Z[l>/V+^(i/ll/l)r"2 

[8]  -‘■O.ZC  2  ]-*-/![  1  2]+ .  *1  ,Bi? 

[9]  'NO  TURNING  POINT .'  ,Z+" 

V 


7  Z  TDP  V \J \N  ;AC\C‘,D 
Cl]  --(((A/lsp^A  1  2  =pp!/)/  3  4), 2 
C 2 ]  -*0 , pl> ' ARGUMENTS  WON "T  DO.' 

C3]  V+(l,f>V)pV 

C4]  ■*•(  ( 0  1  =ppZ)/5,2+4*v/(pZ)='lfplO,2 

[5]  Z-»-(pI/)[2]pZ 

[6]  -*-0-</-*-v/(3=pPPS)a(a/0*  14-11  ,DFS)/^l-ppDPS 

C7]  ->8+(/'/=C/7)A(125>/7)A2SyM+(-/Z?F5C3  1])*0FS[2] 

C  B  J  -0,pO  'DPS  WON ' '  T  DO.' 

C9]  ,DPL=lDPL)*(k=+/pDPL)M-ppDPL 

CIO]  -*11  +  ( l<PPPCl  3  )a>/PPL[  3  l],l>" 

C 11 3  -K),p[ y-' DPI  WON" T  DO.' 

C12]  1  +  ’  ,(N4-\.N)pAC+' . —  +  • 

[13]  Ll:0(/7+2)p'  ' 

C 1 4 ]  C[  1  r  (L  2 . 5+ ( VI  M -DPS[  1  ]  ) tPP5C  2])LV+2>Z[J] 

CIS]  -+Z.3*  1 0 = C 3 ]  |c/-PPL[l] 

C 1 6 ]  •  I  '  ,l+~l+<7 

C17]  h2  :-*£lx  i  (p  IOC  2  ]£</«-</+ 1 

[18]  -H),pl>'  + '  ,NoAC 

[19]  £3  :Z>-(4plO  )tDPL\  2]+PP£[4]*  (</-0PL[  1  ]  )iDPL[  3] 

[20]  -*£2,pl>,01234!il789  '  [1+D+  (  10*(P=0)a  .  v(  1 4  )° .  >,  3 )  ,0] , '  |'  ,U"l+C 
V 
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